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Abstract 



A natural calculus for describing the bound-state structure of relativistic compos- 
ite systems in quantum field theory is the light-front Fock expansion which encodes 
the properties of a hadrons in terms of a set of frame-independent n— particle wave- 
functions. Light-front quantization in the doubly-transverse light-cone gauge has a 
number of remarkable advantages, including explicit unitarity, a physical Fock expan- 
sion, the absence of ghost degrees of freedom, and the decoupling properties needed 
to prove factorization theorems in high momentum transfer inclusive and exclusive 
reactions. A number of applications are discussed in these lectures, including semilep- 
tonic B decays, two-photon exclusive reactions, diffractive dissociation into jets, and 
deeply virtual Compton scattering. The relation of the intrinsic sea to the light-front 
wavefunctions is discussed. Light-front quantization can also be used in the Hamil- 
tonian form to construct an event generator for high energy physics reactions at 
the amplitude level. The light-cone partition function, summed over exponentially- 
weighted light-cone energies, has simple boost properties which may be useful for 
studies in heavy ion collisions. I also review recent work which shows that the struc- 
ture functions measured in deep inelastic lepton scattering are affected by final-state 
rescattering, thus modifying their connection to light-front probability distributions. 
In particular, the shadowing of nuclear structure functions is due to destructive in- 
terference effects from leading-twist diffraction of the virtual photon, physics not 
included in the nuclear light-cone wavefunctions. 



1 Introduction 

Progress in the development and testing of quantum chromodynamics will require 
a detailed understanding of hadron processes at the amplitude level. For example, 
exclusive B-meson decays depend critically on the wavefunction of the B as well as the 
final-state hadronic wavefunctions. Spin correlations such as single-spin asymmetries 
in hard QCD reactions, require an understanding of the phase structure of hadron 
amplitudes, physics well beyond that contained in probability distributions. 

One of the challenges of relativistic quantum field theory is to compute the wave- 
functions of bound states, such as the amplitudes which determine the quark and 
gluon substructure of hadrons in quantum chromodynamics. However, any extension 
of the Heisenberg-Schrodinger formulation of quantum mechanics H\ ip) = i^^\ip) — 
E\'\\}) to the relativistic domain has to confront seemingly intractable hurdles: (1) 
quantum fiuctuations preclude finite particle-number wavefunction representations; 
(2) the charged particles arising from the quantum fiuctuations of the vacuum con- 
tribute to the matrix element of currents - thus knowledge of the wavefunctions alone 
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is insufficient to determine observables; and (3) the boost of an equal-time wavefunc- 
tion from one Lorentz frame to another not only changes particle number, but is as 
complicated a dynamical problem as solving for the wavefunction itself. 

In 1949, Dirac made the remarkable observation that ordinary "instant" time 
t is not the only possible evolution parameter. In fact, evolution in "light-front" time 
T = t + z/c = x~^ has extraordinary advantages for relativistic systems, stemming 
from the fact that 7 out of the 10 Poincare' generators, including a Lorentz boost 
K^, are kinematical (interaction-independent) when one quantizes a theory at fixed 
light-front time. 

The light-front fixes the initial boundary conditions of a composite system as its 
constituents are intercepted by a light- wave evaluated on the hyperplane = t + z/c. 
In contrast, determining an atomic wavefunction at a given instant t = requires 
measuring the simultaneous scattering of Z photons on the Z electrons. In fact, the 
Fock-state representation of bound states defined at equal light-cone time, i.e., along 
the light-front, provides wavefunctions of fixed particle number which are independent 
of the eigenstate's four-momentum P'^. Furthermore, quantum fluctuations of the 
vacuum are absent if one uses light-front time to quantize the system, so that matrix 
elements such as the electromagnetic form factors only depend on the currents of the 
constituents described by the light-cone wavefunctions. I will use here the notation 
Ai" = {A+,A-,A_i_), where A^ = A'^ ± A\ and the metric is A ■ B = liA+B' + 
A-B+) - A^ ■ B^. 

In Dirac's "Front Form", the generator of light-front time translations is P~ = i^- 
Boundary conditions are set on the transverse plane labelled by x± and x~ = z — ct. 
See Fig. |l]. Given the Lagrangian of a quantum field theory, P~ can be constructed 
as an operator on the Fock basis, the eigenstates of the free theory. Since each 
particle in the Fock basis is on its mass shell, k~ = k — k = — , and its energy 
k^ = |(/c+ + A;~) is positive, only particles with positive momenta k'^ = k^ + k^ > can 
occur in the Fock basis. Since the total plus momentum P+ = J2n is conserved, the 
light-cone vacuum cannot have any particle content. The operator H^c = P~^P~ —P]_, 
the "light-cone Hamiltonian" , is frame-independent. 

The Heisenberg equation on the light-front is 

Hlc\^) = M^\^) . (1) 

This can in principle be solved by diagonalizing the matrix (nliJicI^) on the free 
Fock basis: [0] 

{n\HLc\m) (m|V') = {n\^) . (2) 

m 

For example the interaction terms of QCD are illustrated in Fig. The eigenvalues 
{M^} of Hlc = -f^Lc + ^Lc give the squared invariant masses of the bound and 
continuum spectrum of the theory. The light-front Fock space is the eigenstates of 
the free light-front Hamiltonian; i.e., it is a Hilbert space of non-interacting quarks 
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Figure 1: Dirac's three forms of Hamiltonian dynamics. From Ref. 0. 



and gluons, each of which satisfy fc^ = and k~ = "^^^ > 0. The projections 
of the eigensolution on the n-particle Fock states provide the hght-front 
wavef unctions. Thus solving a quantum field theory is equivalent to solving a coupled 
many-body quantum mechanical problem: 

- E -^]i^n = Y.j {n\VLc\n') i^n' (3) 

where the convolution and sum is understood over the Fock number, transverse mo- 
menta, plus momenta, and helicity of the intermediate states. Here m\ = m? + k^. 
An essentially equivalent approach to light-front quantization, pioneered by Weinberg 
P, ^ , is to evaluate the equal-time theory from the perspective of an observer moving 
in the negative z direction with arbitrarily large momentum Pj — oo. The light-cone 
fraction x = ^ of a constituent can be identified with the longitudinal momentum 
X = in a hadron moving with large momentum . Light-front wavefunctions are 
also related to momentum-space Bethe-Salpeter wavefunctions by integrating over 
the relative momenta k~ = — k^ since this projects out the dynamics at = 0. 

We can compare the light-front Fock expansion with the n-particle 
Schrodinger momentum space wavefunction ipN{Pi) of a composite system is the pro- 
jection of the exact eigenstate of the equal-time Hamiltonian on the n-particle states 
of the non-interacting Hamiltonian, the Fock basis. It represents the amplitude for 
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Figure 2: The front-form matrix of QCD interactions in light-cone gauge. Up to 
eight constituents in a meson are shown. From Ref. [Q] and H. C. Pauh. 
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finding the constituents with three- momentum pj, orbital angular momentum, and 
spin, subject to three-momentum conservation and angular momentum sum rules. 
The constituents are on their mass shell, Ei = \Jpi + rnf but do not conserve energy 
Y^^=i Ei> E = ^'p^ + M'^. However, in a relativistic quantum theory, a bound-state 
cannot be represented as a state with a fixed number of constituents. For example, 
the existence of gluons which propagate between the valence quarks necessarily im- 
plies that the hadron wavefunction must describe states with an arbitrary number of 
gluons. Thus a hadronic wavefunction must describe fiuctuations in particle number 
n, as well as momenta and spin. One has to take into account fiuctuations in the 
wavefunction which allow for any number of sea quarks, as long as the total quantum 
numbers of the constituents are compatible with the overall quantum numbers of the 
baryon. 

It is especially convenient to develop the light-front formalism in the light-cone 
gauge = A" + = 0. In this gauge the A~ field becomes a dependent degree 
of freedom, and it can be eliminated from the gauge theory Hamiltonian, with the 
addition of a set of specific instantaneous light-front time interactions. In fact in 
QCD{\ + 1) theory, this instantaneous interaction provides the confining linear x~ 
interaction between quarks. In 3 + 1 dimensions, the transverse field propagates 
massless spin-one gluon quanta with polarization vectors which satisfy both the 
gauge condition = and the Lorentz condition /c ■ e = 0. Thus no extra condition 
on the Hilbert space is required. 

In QCD, the wavefunction of a hadron describes its composition in terms of the 
momenta and spin projections of quark and gluon constituents. For example, the 
eigensolution of a negatively- charged meson QCD, projected on its color-singlet B = 
0, Q = —1, Jz = eigenstates {| ^)} of the free Hamiltonian H^^^{g = 0) at fixed 
T = t — z/c has the expansion: 

n; XiP^, XiP± + k±i, A^J) ipn/Mixi, k±i, Aj). (4) 

The set of light-front Fock state wavefunctions {ipn/Ai} represent the ensemble of 
quark and gluon states possible when the meson is intercepted at the light-front. The 
light-front momentum fractions Xi = kf /P^ = {k^ + k^)/{P^ + P^) with X^ILi^* ~ ^ 
and k_i_i with X^iLi ^-L* ~ Oj- represent the relative momentum coordinates of the QCD 
constituents and are independent of the total momentum of the state. 

Remarkably, the light-front wavefunctions ipn/pi^i, k_i_i, Aj) are independent of the 
proton's momentum P'^ = P^ + P^ , and Pi. Thus once one has solved for the light- 
front wavefunctions, one can compute hadron matrix elements of currents between 
hadronic states of arbitrary momentum. The actual physical transverse momenta 
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are pxi = XiP± + k±i. The Aj label the light-front spin projections of the quarks 
and gluons along the quantization z direction. The spinors of the light-front for- 
malism automatically incorporate the Melosh-Wigner rotation. The physical gluon 
polarization vectors e'^{k, A = ±1) are specified in light-cone gauge by the conditions 
k ■ e = 0, 1] ■ e = e'^ = 0. The parton degrees of freedom are thus all physical; there 
are no ghost or negative metric states. 

The light-front representation thus provides a frame-independent, quantum- mechan- 
ical representation of a hadron at the amplitude level, capable of encoding its multi- 
quark, hidden-color and gluon momentum, helicity, and flavor correlations in the form 
of universal process-independent hadron wavefunctions. 

Angular momentum has simplifying features in the light-front formalism since 
the projection is kinematical and conserved. Each light-front Fock wavefunction 
satisfies the angular momentum sum rule: = J2i=i + Sj=i ^ ■ The sum over 

represents the contribution of the intrinsic spins of the n Fock state constituents. 
The sum over orbital angular momenta 

derives from the n — 1 relative momenta. This excludes the contribution to the orbital 
angular momentum due to the motion of the center of mass, which is not an intrinsic 
property of the hadron. The numerator structure of the light-front wavefunctions is 
in large part determined by the angular momentum constraints. 

If one imposes periodic boundary conditions in x~ = t + z/c, then the plus mo- 
menta become discrete: kf = ^ni,P^ = ^K, where = K [H, 0]. For a 
given "harmonic resolution" K, there are only a finite number of ways positive in- 
tegers Hi can sum to a positive integer K. Thus at a given K, the dimension of 
the resulting light-front Fock state representation of the bound state is rendered fi- 
nite without violating Lorentz invariance. The eigensolutions of a quantum field 
theory, both the bound states and continuum solutions, can then be found by nu- 
merically diagonalizing a frame-independent light-front Hamiltonian Hlc on a finite 
and discrete momentum-space Fock basis. Solving a quantum field theory at fixed 
light-front time r thus can be formulated as a relativistic extension of Heisenberg's 
matrix mechanics. The continuum limit is reached for K —>■ oo. This formulation of 
the non-perturbative light-front quantization problem is called "discretized light-cone 
quantization" (DLCQ) [0]. Lattice gauge theory has also been used to calculate the 
pion light-front wavefunction 

The DLCQ method has been used extensively for solving one-space and one-time 
theories 0, including applications to supersymmetric quantum field theories and 
specific tests of the Maldacena conjecture []TU[. There has been progress in system- 



atically developing the computation and renormalization methods needed to make 
DLCQ viable for QCD in physical spacetime. For example, John Hiller, Gary Mc- 
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Carter, and I |Tl| have shown how DLCQ can be used to solve 3+1 theories despite 
the large numbers of degrees of freedom needed to enumerate the Fock basis. A key 
feature of our work is the introduction of Pauli Villars fields to regulate the UV diver- 
gences and perform renormalization while preserving the frame-independence of the 
theory. A recent application of DLCQ to a 3+1 quantum field theory with Yukawa 
interactions is given in Ref. |Tl[|. There has also been important progress using the 
transverse lattice, essentially a combination of DLCQ in 1+1 dimensions together with 



a lattice in the transverse dimensions O, |T3|, |T^. One can also define a truncated 



theory by eliminating the higher Fock states in favor of an effective potential ||T5 |. 
Spontaneous symmetry breaking and other nonperturbative effects associated with 
the instant-time vacuum are hidden in dynamical or constrained zero modes on the 
light-front. An introduction is given in Refs. |T6|, p!7| . 

Because of their Lorentz invariance, it is particularily easy to write down exact 
expressions for matrix elements of currents and other local operators, even the cou- 
plings of gravitons. In fact as I discuss in Section 3, one can show that the anomalous 
gravito-magnetic moment -B(O), analogous to F2{0) in electromagnetic current in- 
teractions, vanishes identically for any system, composite or elementary [|T^. This 
important feature which follows in general from the equivalence principle, is obeyed 
explicitly in the light-front formalism. 

The set of light-front wavefunctions provide a frame-independent, quantum-mechan- 
ical description of hadrons at the amplitude level capable of encoding multi-quark and 
gluon momentum, helicity, and flavor correlations in the form of universal process- 
independent hadron wavefunctions. Matrix elements of spacelike currents such as 
the spacelike electromagnetic form factors have an exact representation in terms of 
simple overlaps of the light-front wavefunctions in momentum space with the same 
Xi and unchanged parton number |l^, The measurement and interpretation 

of the basic parameters of the electroweak theory and CP violation depends on an 
understanding of the dynamics and phase structure of B decays at the amplitude 
level. The light-front Fock representation is specially advantageous in the study of 
exclusive B decays. For example, we can write down an exact frame-independent 
representation of decay matrix elements such as S — > DiV from the overlap of n' = n 
parton conserving wavefunctions and the overlap oi n' = n — 2 from the annihilation 



of a quark-antiquark pair in the initial wavefunction ||22[. The handbag contribu- 
tion to the leading-twist off-forward parton distributions measured in deeply virtual 
Compton scattering have a similar light-front wavefunction representation as overlap 
integrals of light-front wavefunctions [^, |2^. I will review this application in Sections 
3 and 4. 

Factorization theorems have recently been proven which allow one to rigorously 
compute certain types of exclusive B decays in terms of the light-front wavefunctions 
and distribution amplitudes of B meson and the final state hadrons. The proofs are 
similar to those used in the analysis of exclusive amplitudes involving large momentum 
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transfer. I review this topic in Section 5 and 6. 

In principle, the hght-front wavefunctions contain fluctuations of states with arbi- 
trary number of quark and gluon partons. For example, contains higher Fock states 
such as uudss > and uudcc > which are intrinsic to the physics of the proton itself; 
i.e., they are multi-connected to the valence quarks and are not generated by gluon 
splitting. A rigorous analysis of the momentum fraction and spin carried by intrinsic 



heavy quarks recently been given by Franz et al . These quantities scale nominally 
as I/ttIq in non-Abelian gauge theory, in striking contrast to the ^/ttLq scaling which 
follows from the Euler-Heisenberg Lagrangian in QED. In general, the intrinsic sea in 
the proton is asymmetric between the Q{x) and Q{x) distributions, in contrast to the 
near symmetry of quark and antiquark distributions generated by DGLAP evolution. 

The fact that the B meson contains Fock states with intrinsic strangeness and 
charm leads to a number of new phenomena in exclusive B decays. In particular, 
since the charm quarks can facilitate weak interactions, one can evade the CKM 
hierarchy. Susan Gardner and I have shown that the color octet intrinsic charm Fock 
components of the B meson can give signiflcant modifications of standard predictions 
for channels such as -B ^ pir. I will review this in Section 8. 

The quark and gluon probability distributions qi{x, Q) and g{x, Q) of a hadron can 
be computed from the absolute squares of the light-front wavefunctions, integrated 
over the transverse momentum up to the resolution scale Q. All helicity distributions 
are thus encoded in terms of the light-front wavefunctions. The DGLAP evolution of 
the structure functions can be derived from the high k± properties of the light-front 
wavefunctions. Thus given the light-front wavefunctions, one can compute all of 
the leading twist helicity and transversity distributions measured in polarized deep 
inelastic lepton scattering. For example, the helicity-specific quark distributions at 
resolution A correspond to 



gA,/A.(^, ^) = E / n ^^1^ E k^^, (6) 



n,<?a- j=l J^"" A, 

:lQ7r'6 (l - Ex,) 5(2) (j2 k^) 6{x - x,)6x,xM^' " -^n) 



where the sum is over all quarks which match the quantum numbers, light-front 
momentum fraction x, and helicity of the struck quark. Similarly, the transversity 
distributions and off-diagonal helicity convolutions are defined as a density matrix of 
the light-front wavefunctions. This defines the LC factorization scheme where the 
invariant mass squared = Y17=i (^li + of the n partons of the light-front 

wavefunctions is limited to TVI^ < A^. 

However, it is not true that the leading-twist structure functions Fi{x,Q^) mea- 
sured in deep inelastic lepton scattering are identical to the quark and gluon distri- 
butions. For example, it is usually assumed, following the parton model, that the F2 
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structure function measured in neutral current deep inelastic lepton scattering is at 
leading order in l/Q^ simply F2{x,Q'^) = Z]g Q^), where x = Xbj = Q'^/2p ■ q 

and q{x, Q) can be computed from the absolute square of the proton's light-front 
wavefunction. I will report on recent work by Paul Hoyer, Nils Marchal, Stephane 
Peigne, Francesco Sannino, and myself which shows that this standard identification 
is wrong. In particular, the shadowing corrections related to the Gribov- Glauber 
mechanism, the interference effects of leading twist diffractive processes in nuclei are 
separate effects in deep inelastic scattering, are not computable from the bound state 
wavefunctions of the target nucleon or nucleus. 

Remarkably, it is now possible to measure the light-front wavefunctions of a rela- 
tivistic hadron by diffractively dissociating it into jets whose momentum distribution 
is correlated with the valence quarks' momenta ^ ^ . At high energies each 



light-front Fock state interacts distinctly; e.g., Fock states with small particle num- 
ber and small impact separation have small color dipole moments and can traverse 
a nucleus with minimal interactions. This is the basis for the predictions for "color 
transparency" in hard quasi-exclusive |]30| , ^ and diffractive reactions ^ . 



QGD color transparency thus tests a fundamental ansatz of QCD, that hadronic inter- 
actions are a manifestation of gauge interactions. The E791 experiment has recently 
provided a remarkable confirmation of this consequence of QGD color transparency, a 
key property of LGWFs and the gauge field interactions in QGD. The new EVA spec- 
trometer experiment E850 at Brookhaven has also reported striking effects of color 



transparency in quasi-elastic proton-proton scattering in nuclei . I will review this 
important development in Section 7. 

The GLEG collaboration has verified the scaling and angular predictions for hard 
exclusive two-photon processes such as 7*7 — >■ vr'' and 77 tt+tt^. The L3 experi- 
ment at LEP at GERN has also measured a number of exclusive hadron production 
channels in two-photon processes, providing important constraints on baryon and 
meson distribution amplitudes and checks of perturbative QGD factorization. These 
processes are particularly sensitive to the meson distribution amplitudes, the non- 
perturbative wavefunctions which control hard QGD exclusive processes, information 
essential for progress in interpreting exclusive B decays. New data from GLEO (Paar, 
et al.) for 77 — >■ 7r"'"7r"'" -|- K^K~ sXW = a/s > 2.5 GeV. is in striking agreement with 
the perturbative QGD prediction given by Lepage and myself. Moreover, the angular 
distribution shows a striking transition to the predicted QGD form as W is raised. 
The 7*7 7r° results are in close agreement with the scaling and normalization of 
the PQGD prediction, provided that the pion distribution amplitude Q) is close 
to the x(l — x) form, the asymptotic solution to the evolution equation. In Section 
6 I review the theory and emphasized the need for more such meson pair production 
data, particularly measurements of ratios and angular dependencies which are partic- 
ularly sensitive to the meson and baryon distribution amplitudes 0, (j)M{x,Q), and 
4>B{xi,Q)- These quantities specify how a hadron shares its longitudinal momentum 
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among its valence quarks; they control virtually all exclusive processes involving a 
hard scale Q, including form factors, Compton scattering and photoproduction at 
large momentum transfer, as well as the decay of a heavy hadron into specific final 



states 33, 34 



The discretized light-front quantization method developed by H.C. Pauli and my- 
self 1^51 is a powerful technique for finding the non-perturbative solutions of quantum 
field theories. The basic method is to diagonalize the light-front Hamiltonian in a 
light-front Fock basis defined using periodic boundary conditions in x~ and x±. The 
method preserves the frame-independence of the Front form. The DLCQ method is 
now used extensively to solve one-space and one-time theories, including supersym- 
metric theories. New applications of DLCQ to supersymmetric quantum field theories 
and specific tests of the Maldacena conjecture have recently been given by Pinsky and 
Trittman. 

There has been progress recently in systematically developing the computation 
and renormalization methods needed to make DLCQ viable for QCD in physical 
spacetime. Recently John Hiller, Gary McCartor and I have shown how DLCQ can 
be used to solve 3-1-1 theories despite the large numbers of degrees of freedom needed to 
enumerate the Fock basis [|n[] . A key feature of our work, is the introduction of Pauli 
Villars fields in order to regulate the UV divergences and perform renormalization, 
again while preserving the frame-independence of the theory. Further discussion will 
be given in Section 9. A review of DLCQ and its applications is given in Ref. ||36|| . 
There also has been important progress using the transverse lattice, essentially a 
combination of DLCQ in i-|-l dimensions together with a lattice in the transverse 
space. 

Models of the light-front wavefunction are important in the absence of exact so- 
lutions. A simple but potentially useful model developed by Dae Sung Hwang and 
myself is discussed in Section 10. 

The interaction Hamiltonian of QCD in light-cone gauge can be derived by sys- 
tematically applying the Dirac bracket method to identify the independent fields 
p7| . It contains the usual Dirac interactions between the quarks and gluons, the 



three-point and four-point gluon non-Abelian interactions plus instantaneous light- 
front-time gluon exchange and quark exchange contributions 



TCr, 
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where 

J^a = ^VitaY'r + faUd-A,,)A^'' . (8) 

In light-front time-ordered perturbation theory, a Green's functions is expanded 
as a power series in the interactions with hght-front energy denominators I^initiai k~ — 
S intermediate + replacing the usual energy denominators. [For a review see Ref. 
p8| .] In general each Feynman diagram with n vertices corresponds to the sum of n\ 
time-ordered contributions. However, in light-front-time-ordered perturbation theory, 
only those few graphs where all > survive. In addition the form of the light-front 



kinetic energies is rational: k^ = -^jp: — , replacing the nonanalytic k^ = y k'^ + 
of equal-time theory. Thus light-front-time-ordered perturbation theory provides a 
viable computational method where one can trace the physical evolution of a process. 
The integration measures are only three-dimensional d'^k±dx] in effect, the k~ integral 
of the covariant perturbation theory is performed automatically. 

Alternatively, one derive Feynman rules for QCD in light-cone gauge, thus allow- 
ing the use of standard covariant computational tools and renormalization methods 
including dimensional regularization. Prem Srivastava and I have recently pre- 
sented a systematic study of light-front-quantized gauge theory in light-cone gauge 
using a Dyson- Wick S-matrix expansion based on light-front-time-ordered products. 
The gluon propagator has the form 

(0| nA',(x)A\m |0) = / d*k §^ (9) 

where we have defined 



D^k) = D.,(k) = -g,. + - n,n.. (10) 



Here is a null four-vector, gauge direction, whose components are chosen to be 
= S^^, = 6'^-. Note also 

D^x{k)D\{k)=D^^{k)D^,{k) = -D^, 

k'^D^ = 0, n^D^ = D_,{k) = 0, 

D,,{q)D^''{k)D,,{q') = -D,,{q)D^^^{q'). (11) 

The gauge field propagator i D^^{k)/ {k"^ + ie) is transverse not only to the gauge 
direction but also to fc^, i.e., it is doubly-transverse. This leads to appreciable 
simplifications in the computations in QCD. For example, the coupling of gluons to 
propagators carrying high momenta is automatic. The absence of collinear divergences 
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in irreducible diagrams in the light-cone gauge greatly simplifies the leading-twist 
factorization of soft and hard gluonic corrections in high momentum transfer inclusive 
and exclusive reactions since the numerators associated with the gluon coupling 
only have transverse components. The renormalization factors in the light-cone gauge 
are independent of the reference direction n'^. Since the gluon only has physical 
polarization, its renormalization factors satisfy Zi = Z-^. Because of its explicit 
unitarity in each graph, the doubly-transverse gauge is well suited for calculations 
identifying the "pinch" effective charge |3^, |iO[| . 

The running coupling constant and QCD j3 function have also been computed at 
one loop in the doubly-transverse light-cone gauge |^ . It is also possible to effectively 
quantize QCD using light-front methods in covariant Feynman gauge |^ . 

A remarkable advantage of light-front quantization is that the vacuum state | 0) of 
the full QCD Hamiltonian evidently coincides with the free vacuum. The light-front 
vacuum is effectively trivial if the interaction Hamiltonian applied to the perturbative 
vacuum is zero. Note that all particles in the Hilbert space have positive energy 

= ^{k~^ + k~), and thus positive light-front A;^. Since the plus momenta J2 kf is 
conserved by the interactions, the perturbative vacuum can only couple to states with 
particles in which all kf = 0; i.e., so called zero-mode states. In the case of QED, a 
massive electron cannot have k~^ = unless it also has infinite energy. In a remarkable 
calculation, Bassetto and collaborators |^ have shown that the computation of the 
spectrum of QCD{1 + 1) in equal time quantization requires constructing the full 
spectrum of non perturbative contributions (instantons). In contrast, in the light- 
front quantization of gauge theory, where the k^ = singularity of the instantaneous 
interaction is defined by a simple infrared regularization, one obtains the correct 
spectrum of QCD{1 + 1) without any need for vacuum-related contributions. 

In the case of QCD (3-1-1), the momentum-independent four-gluon non-Abelian 
interaction in principle can couple the perturbative vacuum to a state with four 
collinear gluons in which all of the gluons have all components fcf = 0, thus hinting 
at role for zero modes in theories with massless quanta. In fact, zero modes of auxiliary 
fields are necessary to distinguish the theta-vacua of massless QED(1+1) [^, ^ 
or to represent a theory in the presence of static external boundary conditions or 
other constraints. Zero-modes provide the light-front representation of spontaneous 
symmetry breaking in scalar theories |^ . 



There are other applications of the light-front formalism: 

1. The distribution of spectator particles in the final state in the proton fragmen- 
tation region in deep inelastic scattering at an electron-proton collider are encoded in 
the light-front wavefunctions of the target proton. Conversely, the light-front wave- 
functions can be used to describe the coalescence of comoving quarks into final state 
hadrons. 

2. The light-front wavefunctions also specify the multi-quark and gluon correla- 
tions of the hadron. Despite the many sources of power-law corrections to the deep 
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inelastic cross section, certain types of dynamical contributions will stand out at 
large Xbj since they arise from compact, highly-correlated fluctuations of the proton 
wavefunction. In particular, there are particularly interesting dynamical 
corrections which are due to the interference of quark currents; i.e., contributions 
which involve leptons scattering amplitudes from two different quarks of the target 



nucleon 46 



3. The higher Fock states of the light hadrons describe the sea quark structure 
of the deep inelastic structure functions, including "intrinsic" strangeness and charm 
fluctuations specific to the hadron's structure rather than gluon substructure |l47|, |48[|. 
Ladder relations connecting state of different particle number follow from the QCD 
equation of motion and lead to Regge behavior of the quark and gluon distributions 



at X ^ 49 



4. The gauge- and process-independent meson and baryon valence-quark distribu- 
tion amplitudes (I)m{x,Q), and 4>B{xi,Q) which control exclusive processes involving 
a hard scale Q, including heavy quark decays, are given by the valence light-front 
Fock state wavefunctions integrated over the transverse momentum up to the reso- 
lution scale Q. The evolution equations for distribution amplitudes follow from the 
perturbative high transverse momentum behavior of the light-front wavefunctions 

il- 

5. Proofs of factorization theorems in hard exclusive and inclusive reactions are 
greatly simplified since the propagating gluons in light-cone gauge couple only to 
transverse currents; collinear divergences are thus automatically suppressed. 

6. The deuteron form factor at high is sensitive to wavefunction configura- 
tions where all six quarks overlap within an impact separation b±i < 0{1/Q). The 
leading power-law fall off predicted by QCD is FaiQ"^) = where, 
asymptotically, /(^^(Q^)) oc as{Q'^)^^'^'^' 0, 0. In general, the six-quark wavefunc- 
tion of a deuteron is a mixture of five different color-singlet states. The dominant 
color configuration at large distances corresponds to the usual proton-neutron bound 
state. However at small impact space separation, all five Fock color-singlet compo- 
nents eventually evolve to a state with equal weight, i.e., the deuteron wavefunction 
evolves to 80% "hidden color" . The relatively large normalization of the deuteron 



form factor observed at large hints at sizable hidden-color contributions |^ . Hid- 
den color components can also play a predominant role in the reaction 'yd J/ippn 
at threshold if it is dominated by the multi-fusion process 'jgg ^ J/ip [Q. Hard 
exclusive nuclear processes can also be analyzed in terms of "reduced amplitudes" 
which remove the effects of nucleon substructure. 

Light-front wavefunctions are thus the frame-independent interpolating functions 
between hadron and quark and gluon degrees of freedom. Hadron amplitudes are com- 
puted from the convolution of the light-front wavefunctions with irreducible quark- 
gluon amplitudes. More generally, all multi-quark and gluon correlations in the bound 
state are represented by the light-front wavefunctions. The light-front Fock represen- 
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tation is thus a representation of the underlying quantum field theory. I will discuss 
progress in computing light-front wavefunctions directly from QCD in Sections 9 and 
10. 

Light-front quantization can also be used in the Hamiltonian form to construct an 
event generator for high energy physics reactions at the amplitude level. The light- 
front partition function, summed over exponentially-weighted light-front energies, 
has simple boost properties which may be useful for studies in heavy ion collisions. I 
discuss these topics in Sections 14 and 15. 



2 Other Theoretical Tools 

In addition to the light-front Fock expansion, a number of other useful theoretical 
tools are available to eliminate theoretical ambiguities in QCD predictions: 

(1) Conformal symmetry provides a template for QCD predictions [Q, leading to 
relations between observables which are present even in a theory which is not scale 
invariant. For example, the natural representation of distribution amplitudes is in 
terms of an expansion of orthonormal conformal functions multiplied by anomalous 
dimensions determined by QCD evolution equations |5^, |5^, ^ . Thus an important 
guide in QCD analyses is to identify the underlying conformal relations of QCD 
which are manifest if we drop quark masses and effects due to the running of the 
QCD couplings. In fact, if QCD has an infrared fixed point (vanishing of the Gell- 
Mann-Low function at low momenta), the theory will closely resemble a scale- free 
conformally symmetric theory in many applications. 

(2) Commensurate scale relations ||5^, |5^ are perturbative QCD predictions which 
relate observable to observable at fixed relative scale, such as the "generalized Crewther 



relation" |Q, which connects the Bjorken and Gross-Llewellyn Smith deep inelastic 
scattering sum rules to measurements of the e^e~ annihilation cross section. Such 
relations have no renormalization scale or scheme ambiguity. The coefficients in the 
perturbative series for commensurate scale relations are identical to those of conformal 
QCD; thus no infrared renormalons are present |Q. One can identify the required 
conformal coefficients at any finite order by expanding the coefficients of the usual 
PQCD expansion around a formal infrared fixed point, as in the Banks-Zak method 
^ . All non-conformal effects are absorbed by fixing the ratio of the respective mo- 
mentum transfer and energy scales. In the case of fixed-point theories, commensurate 
scale relations relate both the ratio of couplings and the ratio of scales as the fixed 
point is approached ^ 



(3) ay and Skeleton Schemes. A physically natural scheme for defining the QCD 
coupling in exclusive and other processes is the ayiQ"^) scheme defined from the 
potential of static heavy quarks. Heavy-quark lattice gauge theory can provide highly 
precise values for the coupling. All vacuum polarization corrections due to fermion 
pairs are then automatically and analytically incorporated into the Cell Mann-Low 
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function, thus avoiding the problem of exphcitly computing and resumming quark 
mass corrections related to the running of the coupling ||6ll. The use of a finite 
effective charge such as ay as the expansion parameter also provides a basis for 
regulating the infrared nonperturbative domain of the QCD coupling. A similar 
coupling and scheme can be based on an assumed skeleton expansion of the theory 

mm- 

(4) The Abelian Correspondence Principle. One can consider QCD predictions as 
analytic functions of the number of colors Nc and flavors Np. In particular, one can 
show at all orders of perturbation theory that PQCD predictions reduce to those of 
an Abelian theory at Nc —>■ with a = CpUs and Np = 2Nf/Cf held fixed ||62|| . 
There is thus a deep connection between QCD processes and their corresponding 
QED analogs. 



3 Applications of Light-Front Wavefunctions to 
Current Matrix Elements 

The light-front Fock representation of current matrix elements has a number of simpli- 
fying properties. The space-like local operators for the coupling of photons, gravitons 
and the deep inelastic structure functions can all be expressed as overlaps of light- 
front wavefunctions with the same number of Fock constituents. This is possible since 
one can choose the special frame 
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[] for space-like momentum transfer 
and take matrix elements of "plus" components of currents such as and 
No contributions to the current matrix elements from vacuum fluctuations occur. 
Similarly, given the local operators for the energy- momentum tensor T^'^{x) and the 
angular momentum tensor M^^'^{x), one can directly compute momentum fractions, 
spin properties, and the form factors A{q^) and B{q^) appearing in the coupling of 
gravitons to composite systems 



In the case of a spin-| composite system, the Dirac and Pauli form factors Fi[q^ 
and F2{q^) are defined by 



{P'\rmP) = u{P') \F^{q^)^^ + F2(g 



2M 



u{P) 



(12) 



where q^ = {P' — P)^ and u{P) is the bound state spinor. In the light-front formal- 
ism it is convenient to identify the Dirac and Pauli form factors from the helicity- 
conserving and helicity-flip vector current matrix elements of the current [^] : 



P + qA 



J+(0) 
2P+ 



PA)=F,{q^) 



(13) 



P + qA 



J+(0) 
2P+ 



PA 



iq' 



iq 



2M 



(14) 
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The magnetic moment of a composite system is one of its most basic properties. The 
magnetic moment is defined at the —>■ hmit, 

/^=^[i^i(0)+F2(0)], (15) 

where e is the charge and M is the mass of the composite system. We use the standard 
hght-front frame (g^ = g° ± g^): 

P = {P\P-,P±) = [p\^Ay (16) 

where = —2P ■ q = —q]_ is 4-momentum square transferred by the photon. 

The Pauh form factor and the anomalous magnetic moment n = ^P2{0) can 
then be calculated from the expression 



- (9^ - i^')^ = E / E £*{x., A.) ^i(x., fc^., A.) , (17) 



where the summation is over all contributing Fock states a and struck constituent 
charges ej. The arguments of the final-state light-front wavefunction are 

k±^ = k^i + (1 - Xi)q± (18) 

for the struck constituent and 

k'±i = k±i - Xiq± (19) 

for each spectator. Notice that the magnetic moment must be calculated from the 
spin-flip non-forward matrix element of the current. In the ultra-relativistic limit 
where the radius of the system is small compared to its Compton scale 1/M, the 
anomalous magnetic moment must vanish |^ . The light-front formalism is consistent 
with this theorem. 

The form factors of the energy-momentum tensor for a spin-1 composite are 
defined by 

(P'lT^'^mP) = u(P')[A(g2)7('^r) + i?(g^)^PV)°g. 

+C{q')^{q^q''-9''''q')]u{P) , (20) 

where = (P' - P^, P^ = \{P' + py^ aH"^ = \{a''h'' + a'^fo^), and u{P) is the 
spinor of the system. One can also readily obtain the light-front representation of the 
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A{q ) and B(q ) form factors of the energy-tensor Eq. ( pOD []T8|. In the interaction 
picture, only the non-interacting parts of the energy momentum tensor T++(0) need 
to be computed: 

T++(0) 



P + qA 



r++(o) 



2(P^ 



2M 



(21) 



(22) 



The A(g^) and B{q^) form factors Eqs. (|2l|) and (^) are similar to the Fi{q^) and 
F2{q^) form factors Eqs. (|1^) and ( p!^ with an additional factor of the light-front 
momentum fraction x = / of the struck constituent in the integrand. The B{q^) 
form factor is obtained from the non-forward spin-flip amplitude. The value of -B(O) 
is obtained in the — > limit. The angular momentum projection of a state is given 
by 



yl(0) {U) + [A{Q) + Bmu{P)-a^'u{P) 



(23) 



This result is derived using a wave-packet description of the state. The {U) term 
is the orbital angular momentum of the center of mass motion with respect to an 
arbitrary origin and can be dropped. The coefficient of the {U) term must be 1; 
^4(0) = 1 also follows when we evaluate the four- momentum expectation value {P^)- 
Thus the total intrinsic angular momentum of a nucleon can be identified with the 
values of the form factors A{q^) and -B(g^) at = : 



1 . 

r 



[A(0) + B 



(24) 



The anomalous moment coupling -B(O) to a graviton can in fact be shown to 
vanish for any composite system. This remarkable result, first derived by Okun and 
Kobzarev [63, 34, 65, 66, 57], follows directly from the Lorentz boost properties of 
the light-front Fock representation |]T^ . 

Dae Sung Hwang, Bo-Qiang Ma, Ivan Schmidt, and I [|18| have recently shown that 
the light-front wavefunctions generated by the radiative corrections to the electron in 
QED provides a simple system for understanding the spin and angular momentum 
decomposition of relativistic systems. This perturbative model also illustrates the 
interconnections between Fock states of different number. The model is patterned 
after the quantum structure which occurs in the one- loop Schwinger a/2TX correction 
to the electron magnetic moment In effect, we can represent a spin-| system as 
a composite of a spin-i fermion and spin-one vector boson with arbitrary masses. A 
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Figure 3: Exact representation of electroweak decays and time-like form factors in 
the light-front Fock representation. 



similar model has been used to illustrate the matrix elements and evolution of light- 
front helicity and orbital angular momentum operators [Q. This representation of 
a composite system is particularly useful because it is based on two constituents but 
yet is totally relativistic. We can then explicitly compute the form factors Fi{q^) 
and F2{q^) of the electromagnetic current, and the various contributions to the form 
factors A{q^) and B^q^) of the energy-momentum tensor. 

Another remarkable advantage of the light-front formalism is that exclusive semilep- 
tonic 5-decay amphtudes such as i? ^ A(!m can also be evaluated exactly |2^. The 



time-like decay matrix elements require the computation of the diagonal matrix ele- 
ment n ^ n where parton number is conserved, and the off-diagonal n + 1 — n — 1 
convolution where the current operator annihilates a qq' pair in the initial B wave- 
function. See Fig. |^. This term is a consequence of the fact that the time-like 
decay q^ = {pe + PuY requires a positive light-front momentum fraction g+ > 0. 
Conversely for space-like currents, one can choose = 0, as in the Drell-Yan-West 
representation of the space-like electromagnetic form factors. However, as can be seen 
from the explicit analysis of the form factor in a perturbative model, the off-diagonal 
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convolution can yield a nonzero /q^ limiting form as —>■ 0. This extra term 
appears specifically in the case of "bad" currents such as J~ in which the coupling 
to qq fluctuations in the light-front wavefunctions are favored. In effect, the q~^ —>■ 
limit generates 6{x) contributions as residues of the n + 1 —>■ n — 1 contributions. The 
necessity for such "zero mode" 6{x) terms has been noted by Chang, Root and Yan 
|69| , Burkardt |7^, and Ji and Choi ||7l[| . 

The off-diagonal n + 1 — > n — 1 contributions give a new perspective for the 
physics of -B-decays. A semileptonic decay involves not only matrix elements where 
a quark changes flavor, but also a contribution where the leptonic pair is created 
from the annihilation of a qq' pair within the Fock states of the initial B wavefunc- 
tion. The semileptonic decay thus can occur from the annihilation of a nonvalence 
quark- ant iquark pair in the initial hadron. This feature will carry over to exclusive 
hadronic S-decays, such as — * tt~D~^. In this case the pion can be produced from 
the coalescence of a du pair emerging from the initial higher particle number Fock 
wavefunction of the B. The D meson is then formed from the remaining quarks after 
the internal exchange of a boson. 

In principle, a precise evaluation of the hadronic matrix elements needed for B- 
decays and other exclusive electroweak decay amplitudes requires knowledge of all of 
the light-front Fock wavefunctions of the initial and final state hadrons. In the case 
of model gauge theories such as QCD(1+1) |^ or collinear QCD ||73| in one-space 
and one-time dimensions, the complete evaluation of the light-front wavefunction is 
possible for each baryon or meson bound-state using the DLCQ method. It would be 
interesting to use such solutions as a model for physical i?-decays. 



4 Light-front Representation of Deeply Virtual 
Compton Scattering 

The virtual Compton scattering process ^(7*P 7p) for large initial photon virtu- 
ality = —Q^ has extraordinary sensitivity to fundamental features of the proton's 
structure. Even though the final state photon is on-shell, the deeply virtual process 
probes the elementary quark structure of the proton near the light front as an ef- 
fective local current. In contrast to deep inelastic scattering, which measures only 
the absorptive part of the forward virtual Compton amplitude ImT^*p^^*p, deeply 
virtual Compton scattering allows the measurement of the phase and spin structure 
of proton matrix elements for general momentum transfer squared t. In addition, the 
interference of the virtual Compton amplitude and Bethe-Heitler wide angle scatter- 
ing Bremsstrahlung amplitude where the photon is emitted from the lepton line leads 
to an electron-positron asymmetry in the e^p — > e^7p cross section which is propor- 
tional to the real part of the Compton amplitude fT^, |75|, |76[. The deeply virtual 



Compton amplitude 7*p — > 7p is related by crossing to another important process 
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7*7 ~^ hadron pairs at fixed invariant mass which can be measured in electron-photon 



coUisions fT^] 



To leading order in 1/Q, the deeply virtual Compton scattering amplitude 
7*(g)p(P) l{(l')p{P') factorizes as the convolution in x of the amplitude t^'^ for 
hard Compton scattering on a quark line with the generalized Compton form factors 



H{x, t, C), E{x, t, C), H{x, t, C), and E{x, t, () of the target proton [||, |6|, |78|, ^ |g, 
|8^ , |83| , |8^ , |85| , |86| , p7| . Here x is the light-front momentum fraction of the struck 
quark, and ( = Q'^/2P ■ q plays the role of the Bjorken variable. The square of the 

four-momentum transfer from the proton is given by t = A = 2P-A = —- — 
, where A is the difference of initial and final momenta of the proton (P = P' + A). 
We will be interested in deeply virtual Compton scattering where is large compared 
to the masses and t. Then, to leading order in = ( . Thus ( plays the 

role of the Bjorken variable in deeply virtual Compton scattering. For a fixed value 
of —t, the allowed range of ( is given by 



< c < 



l + i^S - 1 I . (25) 



2M2 i-t) 

The form factor H{x, t, () describes the proton response when the helicity of the pro- 
ton is unchanged, and E{x,t,() is for the case when the proton helicity is flipped. 
Two additional functions H{x,t,Q, and E{x,t,() appear, corresponding to the de- 
pendence of the Compton amplitude on quark helicity. 



Recently, Markus Diehl, Dae Sung Hwang and I |^ have shown how the deeply 
virtual Compton amplitude can be evaluated explicitly in the Fock state represen- 
tation using the matrix elements of the currents and the boost properties of the 
light-front wavefunctions. For the n — n diagonal term (An = 0), the arguments of 
the final-state hadron wavefunction are ^3^, k±i — ^j^A± for the struck quark and 

k±i + ^3^A_L for the n — 1 spectators. We thus obtain formulae for the diagonal 
(parton-number-conserving) contribution to the generalized form factors for deeply 
virtual Compton amplitude in the domain [BH pil ( < Xi < 1: 



c 

T 



I — c 

- C/l(n-.n)(Xl,t,C) - ^y^=^/2(n-^n)(a;i,i,C) 



x<)*(x:, fc^,. A.) V^;^(x„ ku, A.)(yW)^-", (26) 



2^ 



1-C + ^-^l^fnn^n){xM) 
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Figure 4: Light-front time-ordered contributions to deeply virtual Compton scatter- 
ing. Only the contributions of leading twist in are illustrated. These contribu- 
tions illustrate the factorization property of the leading twist amplitude. 
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x<)*(^;, k'^i, A.) A,)(0:^)l-^ (27) 

where 

^1 = 137^ 5 = k±i — i^A_L for the struck quark, 

^ ^ ^ (28) 

= iif^ ) ^li = + i3^^-L for the (n — 1) spectators. 

A sum over all possible hehcities Aj is understood. If quark masses are neglected, the 

— * — * 

currents conserve helicity. Wc also can check that J2^=i = 1) J27=i k'±i — Ox- 

For the n + 1 n — 1 off-diagonal term (An = —2), consider the case where 
partons f and n + 1 of the initial wavefunction annihilate into the current leaving 
n — 1 spectators. Then Xn+i = ( ~ ^i^ ^±n+i = ^± — k±i. The remaining n — 1 
partons have total momentum ((1 — C,)P'^, — Aj_). The final wavefunction then has 
arguments x'^ = and /c^j = k±i + -^A±. We thus obtain the formulae for the 
off-diagonal matrix element of the Compton amplitude in the domain < xi < 



I — 

~ C/l(n+l^«-l)(a:i,t, C) - -^^J==^f2(n+l^n-l){Xi,t,Q 

(Pk^i f (Pk±n+i A /"^ , f (Pk±i 



~ ^Jo 2(27r)3y 2(27r)3 [jjo V 2(27r) 



n+l \ fn+l 



X'0{n-l)K,^li,Ai) '4jl^^^^{{xi,Xi,Xn+l ^ C - Xl}, 

{k±i, k±i, k±n+i — A± — k±i}, {Ai, Aj, A„+i = — Ai}), (29) 



^-^{^^ + 2(1-0 ^ 2M f2{n+l^n-l){xi,t,C) 

^ /-i f (Pk±i_ f(Pk±n+i^ /-i f d^k± i 



n, X"" " "-V-'V " -V-'V i=2 
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{k±i, k±i, k±n+i = A_L - k±i}, {Ai, Ai, A„+i = -Ai}), (30) 

where i = 2,3, ■ ■ ■ ,n label the n — 1 spectator partons which appear in the final-state 
hadron wavefunction with 

= , k'^^ = k^^ + • (31) 

We can again check that the arguments of the final-state wavefunction satisfy Y^^=2 ^'i — 

1, Er=2^i. = o±. 

The above representation is the general form for the generalized form factors of 
the deeply virtual Compton amplitude for any composite system. Thus given the 
light-front Fock state wavefunctions of the eigensolutions of the light-front Hamil- 
tonian, we can compute the amplitude for virtual Compton scattering including all 
spin correlations. The formulae are accurate to leading order in l/Q^. Radiative 
corrections to the quark Compton amplitude of order as{Q'^) from diagrams in which 
a hard gluon interacts between the two photons have also been neglected. 



5 Applications of QCD Factorization to Hard QCD 
Processes 

Factorization theorems for hard exclusive, semi-exclusive, and diffractive processes 
allow the separation of soft non-perturbative dynamics of the bound state hadrons 
from the hard dynamics of a perturbatively-calculable quark-gluon scattering ampli- 
tude. The factorization of inclusive reactions is reviewed in ref. For reviews and 
bibliography of exclusive process calculations in QCD (see Ref. 

The light-front formalism provides a physical factorization scheme which conve- 
niently separates and factorizes soft non-perturbative physics from hard perturbative 
dynamics in both exclusive and inclusive reactions p, pO| . 

In hard inclusive reactions all intermediate states are divided according to A^,^ < 
and Ai"^ > domains. The lower mass regime is associated with the quark and 
gluon distributions defined from the absolute squares of the LC wavefunctions in the 
light front factorization scheme. In the high invariant mass regime, intrinsic transverse 
momenta can be ignored, so that the structure of the process at leading power has the 
form of hard scattering on collinear quark and gluon constituents, as in the parton 
model. The attachment of gluons from the LC wavefunction to a propagator in 
a hard subprocess is power-law suppressed in LC gauge, so that the minimal quark- 
gluon particle-number subprocesses dominate. It is then straightforward to derive the 
DGLAP equations from the evolution of the distributions with logA^. The anomaly 
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contribution to singlet helicity structure function gi{x, Q) can be explicitly identified 
in the LC factorization scheme as due to the '-f*g —>■ qq fusion process. The anomaly 
contribution would be zero if the gluon is on shell. However, if the off-shellness of the 
state is larger than the quark pair mass, one obtains the usual anomaly contribution 
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In exclusive amplitudes, the LC wavefunctions are the interpolating amplitudes 
connecting the quark and gluons to the hadronic states. In an exclusive ampli- 
tude involving a hard scale all intermediate states can be divided according to 
Aif^ < < and < invariant mass domains. The high invariant mass 
contributions to the amplitude has the structure of a hard scattering process Tfj in 
which the hadrons are replaced by their respective (collinear) quarks and gluons. In 
light-cone gauge only the minimal Fock states contribute to the leading power-law 
fall-off of the exclusive amplitude. The wavefunctions in the lower invariant mass 
domain can be integrated up to an arbitrary intermediate invariant mass cutoff A. 
The invariant mass domain beyond this cutoff is included in the hard scattering am- 
plitude Th- The Th satisfy dimensional counting rules |Q. Final-state and initial 



state corrections from gluon attachments to lines connected to the color-singlet dis- 
tribution amplitudes cancel at leading twist. Explicit examples of perturbative QCD 
factorization will be discussed in more detail in the next section. 

The key non-perturbative input for exclusive processes is thus the gauge and 
frame independent hadron distribution amplitude |]90|, |^ defined as the integral of 



the valence (lowest particle number) Fock wavefunction; e.g. for the pion 

M^u A) = / d^k^ ^^/A^u k±^, A) (32) 

where the global cutoff A is identified with the resolution Q. The distribution ampli- 
tude controls leading-twist exclusive amplitudes at high momentum transfer, and it 
can be related to the gauge-invariant Bethe-Salpeter wavefunction at equal light-front 
time. The logarithmic evolution of hadron distribution amplitudes (pH^Xi^Q) can be 
derived from the perturbatively-computable tail of the valence light-front wavefunc- 
tion in the high transverse momentum regime ^. The conformal basis for the 
evolution of the three-quark distribution amplitudes for the baryons ||9^ has recently 



been obtained by V. Braun et al. ||57|| 



The existence of an exact formalism provides a basis for systematic approxima- 
tions and a control over neglected terms. For example, one can analyze exclusive 
semi-leptonic 5-decays which involve hard internal momentum transfer using a per- 
turbative QCD formalism ||9^, |9^, ^ ^ ^ patterned after the perturbative 



analysis of form factors at large momentum transfer. The hard-scattering analysis 
again proceeds by writing each hadronic wavefunction as a sum of soft and hard 
contributions 

= rn\Ml < A') + ^P'^^'^{Ml > A^), (33) 
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where is the invariant mass of the partons in the n-particle Fock state and A is 
the separation scale. The high internal momentum contributions to the wavefunction 
can be calculated systematically from QCD perturbation theory by iterating the 
gluon exchange kernel. The contributions from high momentum transfer exchange to 
the 5-decay amplitude can then be written as a convolution of a hard-scattering 
quark-gluon scattering amplitude with the distribution amplitudes 0(xj,A), the 
valence wavefunctions obtained by integrating the constituent momenta up to the 
separation scale Ain < A < Q. Furthermore in processes such as B ^ ttD where 
the pion is effectively produced as a rapidly-moving small Fock state with a small 
color-dipole interactions, final state interactions are suppressed by color transparency. 



This is the basis for the perturbative hard-scattering analyses |9J, Q ^ 0, . In 
a systematic analysis, one can identify the hard PQCD contribution as well as the 
soft contribution from the convolution of the light-front wavefunctions. Furthermore, 
the hard-scattering contribution can be systematically improved. 

Given the solution for the hadronic wavefunctions tpl^^ with Ai"^ < A^, one can 
construct the wavefunction in the hard regime with Ai^ > A^ using projection opera- 
tor techniques. The construction can be done perturbatively in QCD since only high 
invariant mass, far off-shell matrix elements are involved. One can use this method 
to derive the physical properties of the LC wavefunctions and their matrix elements 

at high invariant mass. Since Ai^ = Y.7=i (^~^'^r~^ ' ^^^^ method also allows the 

derivation of the asymptotic behavior of light-front wavefunctions at large k±, which 
in turn leads to predictions for the fall-off of form factors and other exclusive matrix 
elements at large momentum transfer, such as the quark counting rules for predicting 
the nominal power-law fall-off of two-body scattering amplitudes at fixed 6cm [0 



and helicity selection rules . The phenomenological successes of these rules can be 
understood within QCD if the coupling av{Q) freezes in a range of relatively small 
momentum transfer [p9|l. 



6 Two-Photon Processes 

The simplest and perhaps the most elegant illustration of an exclusive reaction in 
QCD is the evaluation of the photon-to-pion transition form factor 
p.00|| which is measurable in single-tagged two-photon ee eevr^ reactions. The 
form factor is defined via the invariant amplitude F'^ = — ze^F^^(Q^)e^'^^°'p^epgo- • 
As in inclusive reactions, one must specify a factorization scheme which divides the 
integration regions of the loop integrals into hard and soft momenta, compared to the 
resolution scale Q. At leading twist, the transition form factor then factorizes as a 
convolution of the 7*7 — > qq amplitude (where the quarks are collinear with the final 
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Figure 5: (a) Preliminary transverse lattice results for the pion distribution am- 
plitude at _ 10 GeVl The solid curve is the theoretical prediction from the 



combined DLCQ/transverse lattice method ||107|| ; the chain line is the experimental 
result obtained from jet diffractive dissociation ||26[. Both are normalized to the same 
area for comparison, (b) Scaling of the transition photon to pion transition form fac- 
tor Q'^F^^o{Q'^). The dotted and solid theoretical curves are the perturbative QCD 
prediction at leading and next-to-leading order, respectively, assuming the asymptotic 



pion distribution The data are from the CLEO collaboration [|102 



state pion) with the valence light-front wavefunction of the pion: 

F^m(Q') = ^l' dx(l)M{x, Q)T^^_^m{x, Q^). (34) 

The hard scattering amplitude for 77* qq is T!^j^i{x , Q"^) = [{1 — x)Q'^]^^ x 
{1 + 0{as)) . The leading QCD corrections have been computed by Braaten ||101 |. 



The evaluation of the next-to-leading corrections in the physical ay scheme is given 
in Ref. [^. For the asymptotic distribution amplitude 0^'^^™p*(x) = \/3fT^x{l — x) 
one predicts Q^F^^{Q^) = 2/^ (l - |^^^^) where Q* = e'^/^Q is the BLM scale 
for the pion form factor. The PQCD predictions have been tested in measurements 
of 67 eiT^ by the CLEO collaboration ||102|] . See Fig. ^ (b). The observed flat 
scaling of the Q'^F^t,{Q^) data from = 2 to = 8 GeV^ provides an important 
confirmation of the applicability of leading twist QCD to this process. The magnitude 
of Q'^F^t^{Q'^) is remarkably consistent with the predicted form, assuming the asymp- 
totic distribution amplitude and including the LO QCD radiative correction with 
av{e~^^'^Q)/'7T ~ 0.12. One could allow for some broadening of the distribution am- 
plitude with a corresponding increase in the value of ay at small scales. Radyushkin 
||103|| , Ong ||104|| , and Kroll ||105|| have also noted that the scaling and normaliza- 



tion of the photon-to-pion transition form factor tends to favor the asymptotic form 
for the pion distribution amplitude and rules out broader distributions such as the 
two-humped form suggested by QCD sum rules ||lU(j| . 



The two-photon annihilation process 7*7 — > hadrons, which is measurable in 
single-tagged e^e~ — > e^e~hadrons events, provides a semi-local probe of C = -|- 
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hadron systems 7r°, r/", 77', r/c, tt+tt", etc. The 7*7 —>■ tt+tt^ hadron pair process is 
related to virtual Compton scattering on a pion target by crossing. The leading twist 
amplitude is sensitive to the l/x — 1/(1— x) moment of the two-pion distribution 
amplitude coupled to two valence quarks [77| 



Two-photon reactions, 77 HH at large s = (fci + and fixed Ocui, provide a 
particularly important laboratory for testing QCD since these cross-channel "Comp- 
ton" processes are the simplest calculable large-angle exclusive hadronic scattering 
reactions. The helicity structure, and often even the absolute normalization can be 
rigorously computed for each two-photon channel ||100|| . In the case of meson pairs, 
dimensional counting predicts that for large s, s'^da/dt{'y'y MM scales at fixed 
t/s or 6'c.m. up to factors of Ins/A^. The angular dependence of the 77 HH 
amplitudes can be used to determine the shape of the process-independent distribu- 
tion amplitudes, (pniXjQ)- An important feature of the 77 — > MM amplitude for 
meson pairs is that the contributions of Landshoff pitch singularities are power-law 
suppressed at the Born level — even before taking into account Sudakov form factor 
suppression. There are also no anomalous contributions from the x 1 endpoint 
integration region. Thus, as in the calculation of the meson form factors, each fixed- 
angle helicity amplitude can be written to leading order in l/Q in the factorized form 
[Q2 = p2 ^ ^ min(a;Q, {I - x)Q)]: 



M^^^mJj = dx dy(pjj{y, Qy)TH{x, y, s, 9c.m.(pM{x, Q^), (35) 

where is the hard-scattering amplitude 77 — >• {qq){qq) for the production of the 
valence quarks coUinear with each meson, and (I)m{x,Q) is the amplitude for finding 
the valence q and q with light-front fractions of the meson's momentum, integrated 
over transverse momenta k± < Q. The contribution of non- valence Fock states are 
power-law suppressed. Furthermore, the helicity-selection rules [^] of perturbative 
QCD predict that vector mesons are produced with opposite helicities to leading order 
in 1/Q and all orders in a^. The dependence in x and y of several terms in Tx,a' is 
quite similar to that appearing in the meson's electromagnetic form factor. Thus 
much of the dependence on (pAtiXjQ) can be eliminated by expressing it in terms 
of the meson form factor. In fact, the ratio of the 77 — > tt+tt^ and e^e~ fi'^fi~ 
amplitudes at large s and fixed 9cm is nearly insensitive to the running coupling and 
the shape of the pion distribution amplitude: 

f(77^vr+7r-) A\F^{s)\' 

f(77-/^V-) l-COs2^e.m.' ^ ^ 

The comparison of the PQCD prediction for the sum of tt+tt^ plus K^K~ channels 



with recent CLEO data ||108|| is shown in Fig. |g. The CLEO data for charged pion and 
kaon pairs show a clear transition to the scaling and angular distribution predicted by 
PQCD |100| for W = Ms^^ > 2 GeV. See Fig. |. It is clearly important to measure 
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Figure 6: Comparison of the sum of 77 tt+tt^ and 77 K^K~ meson pair 
production cross sections with the scahng and angular distribution of the perturbative 
QCD prediction [[LOCI||. The data are from the CLEO collaboration ||108| . 



the magnitude and angular dependence of the two-photon production of neutral pions 
and p'^p~ cross sections in view of the strong sensitivity of these channels to the shape 
of meson distribution amplitudes. QCD also predicts that the production cross section 
for charged p-pairs (with any helicity) is much larger that for that of neutral p pairs, 
particularly at large 6c.m. angles. Similar predictions are possible for other helicity- 
zero mesons. The cross sections for Compton scattering on protons and the crossed 
reaction 77 — > pp at high momentum transfer have also been evaluated ||109| , |11CI|| , 
providing important tests of the proton distribution amplitude. 

It is particularly compelling to see a transition in angular dependence between 
the low energy chiral and PQCD regimes. The success of leading-twist perturbative 
QCD scaling for exclusive processes at presently experimentally accessible momentum 
transfer can be understood if the effective coupling av{Q*) is approximately constant 
at the relatively small scales Q* relevant to the hard scattering amplitudes The 
evolution of the quark distribution amplitudes In the \ow-Q* domain at also needs to 
be minimal. Sudakov suppression of the endpoint contributions is also strengthened 
if the coupling is frozen because of the exponentiation of a double logarithmic series. 

A debate has continued ||111| , |112| , |113| , |114|| on whether processes such as the pion 
and proton form factors and elastic Compton scattering 7p — > 7p might be dominated 
by higher-twist mechanisms until very large momentum transfer. If one assumes that 
the light-front wavefunction of the pion has the form V'softl^;, k±) = Aexp{—b-^^j^), 
then the Feynman endpoint contribution to the overlap integral at small k± and 
X ~ 1 will dominate the form factor compared to the hard-scattering contribution 
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until very large . However, this ansatz for '?/'soft(a;, has no suppression at fc_L = 
for any x; z.e., the wavefunction in the hadron rest frame does not fall-off at all for 
/c_L = and — oo. Thus such wavefunctions do not represent well soft QCD 

contributions. Endpoint contributions are also suppressed by the QCD Sudakov form 
factor, reflecting the fact that a near-on-shell quark must radiate if it absorbs large 
momentum. One can show |^ that the leading power dependence of the two-particle 
light-front Fock wavefunction in the endpoint region is 1 — x, giving a meson structure 
function which falls as (1 — x)^ and thus by duality a non- leading contribution to the 
meson form factor -F(Q^) cx l/Q^. Thus the dominant contribution to meson form 
factors comes from the hard-scattering regime. 



Radyushkin [112 has argued that the Compton amplitude is dominated by soft 



end-point contributions of the proton wavefunctions where the two photons both 
interact on a quark line carrying nearly all of the proton's momentum. This de- 
scription appears to agree with the Compton data at least at forward angles where 
—t < 10 GeV^. From this viewpoint, the dominance of the factorizable PQCD leading 
twist contributions requires momentum transfers much higher than those currently 
available. However, the endpoint model cannot explain the empirical success of the 
perturbative QCD fixed 6c.m. scaling da / dt{^p tt^u) ~ const at relatively low 
momentum transfer in pion photoproduction [|115|| . 

Clearly much more experimental input on hadron wavefunctions is needed, partic- 
ularly from measurements of two-photon exclusive reactions into meson and baryon 
pairs at the high luminosity B factories. For example, the ratio 
^(77 — > 7r°7r°)/^(77 — > tt+tt") is particularly sensitive to the shape of pion distri- 
bution amplitude. Baryon pair production in two-photon reactions at threshold may 
reveal physics associated with the soliton structure of baryons in QCD | |ll(j| , |117| ]. 
In addition, fixed target experiments can provide much more information on funda- 
mental QCD processes such as deeply virtual Compton scattering and large angle 
Compton scattering. 



7 DifFractive Dissociation and Light-Cone Wave- 
functions 

Diffractive dissociation in QCD can be understood as a three-step process: 

1. The initial hadron can be decomposed in terms of its quark and gluon con- 
stituents in terms of its light-front Fock-state components. 

2. In the second step, the incoming hadron is resolved by Pomeron or Odderon 
(multi-gluon) exchange with the target or by Coulomb dissociation. The exchanged 
interaction has to supply sufficient momentum transfer to put the diffracted state 
X on shell. Light-front energy conservation requires q~ = {nix " ^V)l ^ where mx 
is the invariant mass of X. In a heavy target rest system, the longitudinal momentum 
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transfer is = {fn\ — m^)/i?7riab- Thus the momentum transfer t = to the target 
can be sufficiently small so that the target remains intact. 

In perturbative QCD, the pomeron is generally be represented as multiple gluon 
exchange between the target and projectile. Effectively this interaction occurs over 
a short light-front time interval, and thus like photon exchange, the perturbative 
QCD pomeron can be effectively represented as a local operator. This description 
is believed to be applicable when the pomeron has to resolve compact states and is 
the basis for the terminology "hard pomeron" . The BFKL formalism generalizes the 
perturbative QCD treatment by an all-orders perturbative resummation, generating a 
pomeron with a fixed Regge intercept ap(0). Next to leading order calculations with 
BLM scale fixing leads to a predicted intercept ap(0) ~ 0.4 [|118|| . However, when 



the exchange interactions are soft, a multiperipheral description in terms of meson 
ladders may dominate the physics. This is the basis for the two-component pomeron 
model of Donnachie and Landshoff ||119 . 



Consider a coUinear frame where the incident momentum is large and s = 
{p-K + Ptarget)^ — ptPta.rget- The matrix element of an exchanged gluon with momen- 
tum Qi between the projectile and an intermediate state | A^) is dominated by the 
"plus current": ^7r|j+(0) exp{i^q^ x~ — iq±i ■ x±\nY Note that the coherent sum of 
couplings of an exchanged gluon to the pion system vanishes when its momentum is 
small compared to the characteristic momentum scales in the projectile light-front 
wavefunction: q-^'-Ax^ <^ 1 and qfAx~ -C 1. The destructive interference of the 
gauge couplings to the constituents of the projectile follows simply from the fact that 
the color charge operator has zero matrix element between distinct eigenstates of the 
QCD Hamiltonian: {A\Q\B) = J £x±dx- {A\j+{0)\B) = |T2g. At high energies 
the change in k^' of the constituents can be ignored, so that Fock states of a hadron 
with small transverse size interact weakly even in a nuclear target because of their 



small dipole moment |^0|, ^7j. To a good approximation the sum of couplings to 
the constituents of the projectile can be represented as a derivative with respect to 
transverse momentum. Thus photon exchange measures a weighted sum of transverse 
derivatives dk^ipni^i, k±^, Aj), and two-gluon exchange measures the second transverse 
partial derivative ||121 . 



3. The final step is the hadronization of the n constituents of the projectile Fock 
state into final state hadrons. Since ql' is small, the number of partons in the initial 
Fock state and the final state hadrons are unchanged. Their coalescence is thus 
governed by the convolution of initial and final-state Fock state wavefunctions. In 
the case of states with high k±, the final state will hadronize into jets, each refiecting 
the respective Xi of the Fock state constituents. In the case of higher Fock states 
with intrinsic sea quarks such as an extra cc pair (intrinsic charm), one will observe 
leading J/ip or open charm hadrons in the projectile fragmentation region; i.e., the 
hadron's fragments will tend to have the same rapidity as that of the projectile. 

For example, diffractive multi-jet production in heavy nuclei provides a novel way 
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to measure the shape of the LC Fock state wavefunctions and test color transparency. 
Consider the reaction ^ |122|| ttA —>■ Jeti + Jet2 + A' at high energy where the 
nucleus A' is left intact in its ground state. The transverse momenta of the jets 
balance so that k±i + k±2 = < R~^a ■ The light-front longitudinal momentum 
fractions also need to add to xi + X2 ~ 1 so that ApL < R]^ ■ The process can 
then occur coherently in the nucleus. Because of color transparency, the valence 
wavefunction of the pion with small impact separation, will penetrate the nucleus 
with minimal interactions, diffracting into jet pairs p^- The Xi = x, = 1 



X 



dependence of the di-jet distributions will thus reflect the shape of the pion valence 
light-front wavefunction in x; similarly, the fc_Li — fc_L2 relative transverse momenta 
of the jets gives key information on the second derivative of the underlying shape 
of the valence pion wavefunction I^B], |122| , |121|| . The diffractive nuclear amplitude 



extrapolated to t = should be linear in nuclear number A if color transparency is 
correct. The integrated diffractive rate should then scale as A? / B\ ~ A^/^. 

The results of a diffractive dijet dissociation experiment of this type E791 at 
Fermilab using 500 GeV incident pions on nuclear targets ||123|| appear to be consistent 
with color transparency. The measured longitudinal momentum distribution of the 
jets ||124|| is consistent with a pion light-front wavefunction of the pion with the shape 



of the asymptotic distribution amplitude, 0^^y™P*( 



X 



v^/^x(l — x). Data from 



CLEO pi for the 77* 



transition form factor also favor a form for the pion 



distribution amplitude close to the asymptotic solution to the perturbative QCD 
evolution equation p. 

The interpretation of the diffractive dijet processes as measures of the hadron 
distribution amplitudes has recently been questioned by Braun et al. |125| and by 



Chernyak ||126|| who have calculated the hard scattering amplitude for such processes 
at next-to-leading order. However, these analyses neglect the integration over the 
transverse momentum of the valence quarks and thus miss the logarithmic ordering 
which is required for factorization of the distribution amplitude and color-filtering in 
nuclear targets. 

As noted above, the diffractive dissociation of a hadron or nucleus can also oc- 
cur via the Coulomb dissociation of a beam particle on an electron beam {e.g. at 
HERA or eRHIC) or on the strong Coulomb field of a heavy nucleus {e.g. at RHIC 
or nuclear collisions at the LHC) ||121 ]. The amplitude for Coulomb exchange at 
small momentum transfer is proportional to the first derivative of the light- 

front wavefunction, summed over the charged constituents. The Coulomb exchange 
reactions fall off less fast at high transverse momentum compared to pomeron ex- 
change reactions since the light-front wavefunction is effective differentiated twice in 
two-gluon exchange reactions. 

It will also be interesting to study diffractive tri-jet production using proton beams 
pA Jeti + Jet2 + Jets + A' to determine the fundamental shape of the 3-quark struc- 
ture of the valence light-front wavefunction of the nucleon at small transverse sepa- 
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ration |2^. For example, consider the Coulomb dissociation of a high energy proton 
at HERA. The proton can dissociate into three jets corresponding to the three-quark 
structure of the valence light-front wavefunction. We can demand that the produced 
hadrons all fall outside an opening angle 9 in the proton's fragmentation region. Ef- 
fectively all of the light-front momentum Yl,j — 1 of the proton's fragments will 
thus be produced outside an "exclusion cone" . This then limits the invariant mass of 
the contributing Fock state A^^ > = P^^ sin^ 9/4 from below, so that perturbative 
QCD counting rules can predict the fall-off in the jet system invariant mass Ai. At 
large invariant mass one expects the three-quark valence Fock state of the proton to 
dominate. The segmentation of the forward detector in azimuthal angle (p can be 
used to identify structure and correlations associated with the three-quark light-front 



wavefunction ||121|| . An interesting possibility is that the distribution amplitude of 



the A(1232) for Jz = 1/2,3/2 is close to the asymptotic form X1X2X3, but that the 
proton distribution amplitude is more complex. This ansatz can also be motivated 
by assuming a quark-diquark structure of the baryon wavef unctions. The differences 
in shapes of the distribution amplitudes could explain why the p A transition 
form factor appears to fall faster at large than the elastic p ^ p and the other 



p N* transition form factors ||127|| . One can use also measure the dijet structure 
of real and virtual photons beams 'j*A — >■ Jeti -|- Jet2 -l- A' to measure the shape 
of the light-front wavefunction for transversely-polarized and longitudinally-polarized 
virtual photons. Such experiments will open up a direct window on the amplitude 
structure of hadrons at short distances. The light-front formalism is also applicable 
to the description of nuclei in terms of their nucleonic and mesonic degrees of free- 
dom [[128| , |129|| . Self- resolving diffractive jet reactions in high energy electron-nucleus 
collisions and hadron-nucleus collisions at moderate momentum transfers can thus be 
used to resolve the light-front wavefunctions of nuclei. 

The first tests of color transparency involved large momentum transfer quasi- 
elastic scattering processes in nuclei. Such reactions are predicted in perturbative 
QCD to depend on the scattering of small impact size hadron wavefunction configu- 
rations |^D|]. The onset of color transparency in proton-proton scattering in nuclei was 
first seen by Experiment E834 at BNL by observing a rise in the ratio of quasi-elastic 
to elastic pp scattering at large angles and energies up to i/i ~ 5 GeV | ]13(J| | . Quasi- 
elastic proton-proton scattering is advantageous over the analogous electron-proton 
scattering reaction since the wavefunctions of the incoming and outgoing hadron in 
high energy proton reactions would not suffer rapid expansion. However, E834 also 
revealed another remarkable feature of quasi-elastic pp scattering: the quenching of 
color transparency at the largest measured energy measured by E834, in direct con- 
tradiction to the predictions perturbative QCD. A more recent experiment, E850, 
using the EVA spectrometer has now confirmed this unexpected effect through new 
measurements of the transparency ratio at higher energies [^ . 

The quenching of color transparency observed in the E834 and E850 experiments 
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is almost as important discovery as color transparency itself. It signals a nonpertur- 
bative effect in QCD which clearly must be understood. The quenching occurs at 
the center-of-mass energy of 5 GeV where the pp elastic cross section also displays 
another remarkable effect: the rate of scattering where the spins of the initial protons 
are parallel and normal to the scattering plane grows rapidly and becomes approxi- 
mately 4 times as large as the spin-antiparallel rate ||131|| . De Teramond and I [|132 



have noted that both phenomenon occur just at the threshold for open charm hadron 
production. We have shown in fact that resonance production in pp elastic scattering 
due to a uuduudcc spin-1 resonance will in fact lead to a remarkably large spin cor- 
relation Atvat and quenching of color transparency above the charm threshold. If this 
explanation is validated (by the observation of a significant open charm cross section 
near 5 GeV center of mass energy), then E834 and E850 will have provided the first 
evidence for an exotic QCD state with hidden charm. 



8 Higher Fock States and the Intrinsic Sea 

Since a hadronic wavefunction describes states off of the light-front energy shell, there 
is a finite probability of the projectile having fluctuations containing extra quark- 
antiquark pairs, such as intrinsic strangeness charm, and bottom. In contrast to the 
quark pairs arising from gluon splitting, intrinsic quarks are multiply-connected to 
the valence quarks and are thus part of the dynamics of the hadron. 

Recently Franz, Polyakov, and Goeke have analyzed the properties of the intrinsic 
heavy-quark fluctuations in hadrons using the operator-product expansion |^5[. For 
example, the light-cone momentum fraction carried by intrinsic heavy quarks in the 
proton Xqq as measured by the T~^~^ component of the energy-momentum tensor is re- 
lated in the heavy-quark limit to the forward matrix element {p\tTc{G^°'G'^^Gaf3) /ttlqIp) , 
where C^" is the gauge field strength tensor. Diagrammatically, this can be described 
as a heavy quark loop in the proton self-energy with four gluons attached to the 
light, valence quarks. Since the non-Abelian commutator [Aq,,A^] is involved, the 
heavy quark pairs in the proton wavefunction are necessarily in a color-octet state. 
It follows from dimensional analysis that the momentum fraction carried by the QQ 
pair scales as k^/uiQ where k±_ is the typical momentum in the hadron wave function. 
[In contrast, in the case of Abelian theories, the contribution of an intrinsic, heavy 
lepton pair to the bound state's structure first appears in 0{l/m\). One relevant op- 
erator corresponds to the Born-Infeld {F^^Y light-by- light scattering insertion, and 
the momentum fraction of heavy leptons in an atom scales as fc^/m|^.] 

The intrinsic sea is thus sensitive to the hadronic bound-state structure ||133| , |47[] . 



The maximal contribution of an intrinsic heavy quark occurs at xq ~ rnxq/ J2i''TT'± 



where m± = ^Jm'^ + k'j_] i.e. at large xq, since this minimizes the invariant mass 
A4n- The measurements of the charm structure function by the EMC experiment 
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are consistent with intrinsic charm at large x in the nucleon with a probabihty of 
order 0.6 ±0.3% which is consistent with the recent estimates based on instanton 
fluctuations [p5[] . 

Thus one can identify two contributions to the heavy quark sea, the "extrinsic" 
contributions which correspond to ordinary gluon splitting, and the "intrinsic" sea 
which is multi-connected via gluons to the valence quarks. Intrinsic charm can be 
materialized by diffractive dissociation into open or hidden charm states such as 
pp — > J /ipXp' , KcXp' . At HERA one can measure intrinsic charm in the proton by 
Coulomb dissociation: pe —>■ AcXe', and J/ipXe'. Since the intrinsic heavy quarks 
tend to have the same rapidity as that of the projectile, they are produced at large 
Xi? in the beam fragmentation region. 

The presence of intrinsic charm quarks in the B wave function provides new 
mechanisms for B decays. The characteristic momenta characterizing the B meson is 
most likely higher by a factor of 2 compared to the momentum scale of light mesons. 
This effect is analogous to the higher momentum scale of muonium fi~^e~ versus that 
of positronium e"'"e~ in atomic physics because of the larger reduced mass. Thus one 
can expect a higher probability for intrinsic charm in heavy hadrons compared to 
light hadrons. For example, Chang and Hou have considered the production of final 
states with three charmed quarks such as S — >■ J/ipDir and B J/ipD* [ |134] ] ; these 
final states are difficult to realize in the valence model, yet they occur naturally when 
the b quark of the intrinsic charm Fock state | bucc) decays via b cud. In fact, 
the J/tp spectrum for inclusive B J/tpX decays measured by CLEO and Belle 
shows a distinct enhancement at the low J/ip momentum where such decays would 
kinematically occur. Alternatively, this excess could reflect the opening of baryonic 
channels such as B ^ J/ippA [p.35|] . 

Recently, Susan Gardner and I have shown that the presence of intrinsic charm in 
the hadrons' light-front wave functions, even at a few percent level, provides new, com- 
petitive decay mechanisms for B decays which are nominally CKM-suppressed | |136| ]. 
For example, the weak decays of the S-meson to two-body exclusive states consisting 
of strange plus light hadrons, such as B —y ttK, are expected to be dominated by 
penguin contributions since the tree-level b —>■ suu decay is CKM suppressed. How- 
ever, higher Fock states in the B wave function containing charm quark pairs can 
mediate the decay via a CKM-favored b sec tree-level transition. Such intrinsic 
charm contributions can be phenomenologically significant. Since they mimic the am- 



plitude structure of "charming" penguin contributions ||137|, |138|, |139|, |140||, charming 



penguins need not be penguins at all ||13(: 

One can also distinguish "intrinsic gluons" [|141|| which are associated with multi- 
quark interactions and extrinsic gluon contributions associated with quark substruc- 
ture. One can also use this framework to isolate the physics of the anomaly contri- 



bution to the Ellis- Jaffe sum rule |9T|]. Thus neither gluons nor sea quarks are solely 



generated by DGLAP evolution, and one cannot define a resolution scale Qq where 
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the sea or gluon degrees of freedom can be neglected. 

It is usually assumed that a heavy quarkonium state such as the J/V' always 
decays to light hadrons via the annihilation of its heavy quark constituents to gluons. 
However, as Karliner and I ||142|| have shown, the transition J /ip ^ can also occur 
by the rearrangement of the cc from the J /ip into the | qqcc) intrinsic charm Fock state 
of the p or TT. On the other hand, the overlap rearrangement integral in the decay 
ip' piT will be suppressed since the intrinsic charm Fock state radial wavefunction 
of the light hadrons will evidently not have nodes in its radial wavefunction. This 
observation provides a natural explanation of the long-standing puzzle ||143|| why 
the J/ip decays prominently to two-body pseudoscalar- vector final states, breaking 
hadron helicity conservation whereas the ip' does not. 

The higher Fock state of the proton | uudss) should resemble a | KK) intermediate 
state, since this minimizes its invariant mass M.. In such a state, the strange quark 
has a higher mean momentum fraction x than the s ||144| , |145| , |146|| . Similarly, the 
helicity of the intrinsic strange quark in this configuration will be anti-aligned with 
the helicity of the nucleon | |144| , |146|| . This Q ^ Q asymmetry is a striking feature of 
the intrinsic heavy-quark sea. 



9 Non-Perturbative Solutions of Light-Front Quan- 
tized QCD 

Is there any hope of computing light-front wavefunctions from first principles? The 
solution of the light-front Hamiltonian equation H^^^\ \E') = Af^| \I') is an eigenvalue 
problem which in principle determines the masses squared of the entire bound and 
continuum spectrum of QCD. If one introduces periodic or anti-periodic boundary 
conditions, the eigenvalue problem is reduced to the diagonalization of a discrete 
Hermitian matrix representation of H^q^ ■ The light-front momenta satisfy 
and = ^K, where J^i^i = K. The number of quanta in the contributing Fock 
states is restricted by the choice of harmonic resolution. A cutoff on the invariant mass 
of the Fock states truncates the size of the matrix representation in the transverse 



momenta. This is the essence of the DLCQ method ||35|, which has now become a 
standard tool for solving both the spectrum and light-front wavefunctions of one-space 
one-time theories — virtually any 1-1-1 quantum field theory, including "reduced QCD" 
(which has both quark and gluonic degrees of freedom) can be completely solved 
using DLCQ [|147| , [73| . The method yields not only the bound-state and continuum 



spectrum, but also the light-front wavefunction for each eigensolution ||148| , |149| . 

Dalley et al. have shown how one can use DLCQ in one space-one time, with 
a transverse lattice to solve mesonic and gluonic states in 3 + 1 QCD ||150||. The 



spectrum obtained for gluonium states is in remarkable agreement with lattice gauge 
theory results, but with a huge reduction of numerical effort. Hiller and I [|151|] have 
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shown how one can use DLCQ to compute the electron magnetic moment in QED 
without resort to perturbation theory. 

There has been recent progress developing the computational tools and renormal- 
ization methods which can make DLCQ a viable computational method for QCD in 
physical space-time. John Hiller, Gary McCartor, and I |11| have shown how DLCQ 
can be used to solve (3+1) theories and obtain the spectrum and light-front wave- 
functions of the bound state solutions despite the large numbers of degrees of freedom 
needed to enumerate the Fock basis. A key feature of our work is the introduction 



of Pauli Villars fields [[152|] in the DLCQ basis which regulate the UV divergences 
and perform renormalization while preserving the frame-independence of the theory 

Bin- 



A recent application of DLCQ and Pauli Villars regularization to a (3+1) quantum 
field theory with Yukawa interactions is given in Ref. [11|. Only one heavy fermion 
is allowed in the Fock states. We include an additional effective interaction which 
represents the contribution of the missing Z graph and cancels an infrared singularity 
introduced by the instantaneous fermion interaction. Cancellation of ultraviolet in- 
finities is then arranged by choosing imaginary couplings or an indefinite metric. In 
our most recent work we used three heavy scalars, two of which have negative norm. 

In DLCQ, all light-front momentum variables are discretized, with p+ —>■ nn/L 
and p± n_i_7T/L±, in terms of longitudinal and transverse length scales L and L±. 
The total longitudinal momentum is = Kn/L and momentum fractions are given 
hj X = n/K. Wave functions and the mass eigenvalue problem, where i^Lc = P^P~ ■, 
are naturally expressed in terms of momentum fractions and the resolution K. Hence 
L disappears, and K effectively takes its place as the resolution scale. The transverse 
scale L_[_ is set by a momentum cutoff and a transverse resolution. The integrals over 
wave functions which make up the mass eigenvalue problem Hi^c^ = are then 

approximated by the trapezoidal quadrature rule. This yields a matrix eigenvalue 
problem which is typically quite large but also quite sparse. Lanczos techniques [ |155|| 
are used to extract eigenvalues and eigenvectors for the lowest states, even in the case 
of an indefinite metric \T1 . 

The mass M of the dressed single-fermion state is held fixed. This is imposed by 
rearranging the mass eigenvalue problem into an eigenvalue problem for the quantity 



X 



2 M^+pI ^fi^ + ql, 



X 



(37) 



where /C represents the original kernel and amplitudes are related by </) = ^/xcf). 

The coupling g is constrained by imposing a condition on the boson occupation 
number: (:0^(O):) = $J. :0^(O): $o-- This quantity can be computed fairly efficiently as 
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Figure 7: The boson structure function fs at various numerical resolutions for 
(:0^(O):) = 0.5, with M = fi, cutoff = and Pauli-Villars masses /if = 

20yU^, and /i| = 30fi^. The solid line is from first-order perturbation theory. 



the sum 



CO „ ntot 



(3^ 



71i=0 j 



The constraint on (:0^(O):) can be satisfied by solving it simultaneously with the 
eigenvalue problem. 

With the parameters fixed, we can compute various quantities, such as the parton 
wavefunctions and momentum distributions, the form factor slope at zero momentum 
transfer, the average numbers of constituents, and the average constituent momenta. 
A representative plot of the bosonic structure function 



^, . ^ 7 -1 



ni=0 s 

x%-g+/^+) 



k=l 
2 



(39) 



is given in Fig. |^. We have also worked at somewhat stronger couplings where devia- 
tions from first order perturbation theory becomes apparent; however, high resolution 
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is required, with = 21 to 39 and as many as 15 transverse momentum points. This 
resolution could be achieved by limiting the number of constituents to 3, after veri- 
fying that the contribution from higher sectors was sufficiently small. 

The success of application of DLCQ to the Yukawa theory with Pauli-Villars regu- 
larization is encouraging. One can compute masses and wave functions for eigenstates 
for quantum field theories in physical space-time. Another procedure, now under in- 
vestigation, is to use one heavy scalar and one heavy fermion, both with negative 
norm, as suggested by the work of Paston et al. [|156|| This method has the advantage 
of being free of instantaneous fermion interactions. An alternative and interesting 
regularization is to apply DLCQ to finite supersymmetric 3+1 theories, and then 
introduce supersymmetric breaking. 

We plan to continue to be explored these possibilities with various model theories, 
leading eventually to the direct application to QCD(3+1). In fact, Paston et al. [|157 



have already obtained a PV-like regularization of QCD which could, in principle, be 
solved by DLCQ; however, given present computing power the number of fields is 
possibly too large for meaningful calculations. 

One can also formulate DLCQ so that supersymmetry is exactly preserved in the 
discrete approximation, thus combining the power of DLCQ with the beauty of su- 
persymmetry ||158| , |159| , 116011 . The "SDLCQ" method has been applied to several 



interesting supersymmetric theories, to the analysis of zero modes, vacuum degen- 
eracy, massless states, mass gaps, and theories in higher dimensions, and even tests 



of the Maldacena conjecture ||158|| . Broken supersymmetry is interesting in DLCQ, 



since it may serve as a method for regulating non-Abelian theories [p. 54] . 

There are also many possibilities for obtaining approximate solutions of light- 
front wavefunctions in QCD. QCD sum rules, lattice gauge theory moments, and QCD 
inspired models such as the bag model, chiral theories, provide important constraints. 
Guides to the exact behavior of LC wavefunctions in QCD can also be obtained from 
analytic or DLCQ solutions to toy models such as "reduced" QCD{1 + 1). The light- 
front and many-body Schrodinger theory formalisms must match In the nonrelativistic 
limit. 

It would be interesting to see if light-front wavefunctions can incorporate chiral 
constraints such as soliton (Skyrmion) behavior for baryons and other consequences 
of the chiral limit in the soft momentum regime. Solvable theories such as QCD{1 + 
1) are also useful for understanding such phenomena. It has been shown that the 
anomaly contribution for the 7r° — > 77 decay amplitude is satisfied by the light-front 
Fock formalism in the limit where the mass of the pion is light compared to its size 
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10 Non-Perturbative Calculations of the Pion Dis- 
tribution Amplitude 

The distribution amplitude Q) can be computed from the integral over transverse 
momenta of the renormalized hadron valence wavefunction in the light-cone gauge at 
fixed light-front time ||38|| : 



Q) = J d^kl 9 I^Q' - ^^j^^ i^^'^Xx, k[), (40) 

where a global cutoff in invariant mass is identified with the resolution Q. The distri- 
bution amplitude 0(x, Q) is boost and gauge invariant and evolves in InQ through an 
evolution equation ^. Since it is formed from the same product of operators 

as the non-singlet structure function, the anomalous dimensions controlling 0(x, Q) 
dependence in the ultraviolet logQ scale are the same as those which appear in the 
DGLAP evolution of structure functions |5^. The decay tt — > /xi/ normalizes the 
wave function at the origin: ao/6 = Jq dx(f){x, Q) = f.,r/{2y/3). One can also compute 
the distribution amplitude from the gauge invariant Bethe-Salpeter wavefunction at 
equal light-front time. This also allows contact with both QCD sum rules and lattice 
gauge theory; for example, moments of the pion distribution amplitudes have been 
computed in lattice gauge theory ||162| , |163| , |164| . 



Dalley ||107|| has recently calculated the pion distribution amplitude from QCD 
using a combination of the discretized DLCQ method for the x~ and light-front co- 
ordinates with the transverse lattice method ||165| , |166|| in the transverse directions, A 



finite lattice spacing a can be used by choosing the parameters of the effective theory 
in a region of renormalization group stability to respect the required gauge, Poincare, 
chiral, and continuum symmetries. The overall normalization gives = 101 MeV 
compared with the experimental value of 93 MeV. Figure |^ (a) compares the result- 
ing DLCQ/transverse lattice pion wavefunction with the best fit to the diffractive 
di-jet data (see the next section) after corrections for hadronization and experimental 
acceptance p6|. The theoretical curve is somewhat broader than the experimental 
result. However, there are experimental uncertainties from hadronization and theo- 
retical errors introduced from finite DLCQ resolution, using a nearly massless pion, 
ambiguities in setting the factorization scale Q^, as well as errors in the evolution 
of the distribution amplitude from 1 to 10 GeV^. Instanton models also predict a 



pion distribution amplitude close to the asymptotic form ||167|| . In contrast, recent 



lattice results from Del Debbio et al. ||164|| predict a much narrower shape for the 



pion distribution amplitude than the distribution predicted by the transverse lattice. 
A new result for the proton distribution amplitude treating nucleons as chiral solitons 



has recently been derived by Diakonov and Petrov [168|. Dyson- Schwinger models 



[|169|] of hadronic Bethe-Salpeter wavefunctions can also be used to predict light-front 
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wavef unctions and hadron distribution amplitudes by integrating over the relative 
k~ momentum. There is also the possibility of deriving Bethe-Salpeter wavefunc- 



tions within light-front gauge quantized QCD [|T| in order to properly match to the 
light-cone gauge Fock state decomposition. 

11 Calculating and Modelling Light-Front Wave- 
functions 

Many features of the light-front wavefunctions follow from general arguments. Light- 
front wavefunctions satisfy the equation of motion: 

H^^^'l^) = {Hl^ + V,c)\^) = M'\^), 

which has the Heisenberg matrix form in Fock space: 



n 2 

[M'-j:^]ij^ = Y.[{n\V\n')A, 

1=1 n' 



where the convolution and sum is understood over the Fock number, transverse mo- 
menta, plus momenta and helicity of the intermediate states. Here m\ = im? + k^. 
Thus, in general, every light-front Fock wavefunction has the form: 



n 



2 

±2 



where r„ = J2n' I ynn'^^n- The main dynamical dependence of a light-front wave- 
function away from the extrema is controlled by its light-front energy denominator. 
The maximum of the wavefunction occurs when the invariant mass of the partons is 
minimal; i.e., when all particles have equal rapidity and are all at rest in the rest 
frame. In fact, Dae Sung Hwang and I ||12CI|| have noted that one can rewrite the 
wavefunction in the form: 

/ r„ 

tPn 



where Xi = Xi = mxi/ Yl^=ifn±_i is the condition for minimal rapidity differences of 
the constituents. The key parameter is — Y^^=i^\i/ = — M^(5^. We can also 
interpret 5^ ~ 2e/M where e = Z^ILi "^J-j — ^ is the effective binding energy. This 
form shows that the wavefunction is a quadratic form around its maximum, and that 
the width of the distribution in {xi — XiY (where the wavefunction falls to half of 
its maximum) is controlled by and the transverse momenta Note also that 
the heaviest particles tend to have the largest Xj, and thus the largest momentum 
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fraction of the particles in the Fock state, a feature famihar from the intrinsic charm 
model. For example, the b quark has the largest momentum fraction at small k_i_ in 
the B meson's valence light-front wavefunction,, but the distribution spreads out to 
an asymptotically symmetric distribution around Xb ~ 1/2 when k± 3> ml. 

We can also discern some general properties of the numerator of the light-front 
wavef unctions. r„(xj, fc_Lj, Aj). The transverse momentum dependence of r„ guar- 
antees Jz conservation for each Fock state. For example, one of the three light- 
front Fock wavefunctions of a = +1/2 lepton in QED perturbation theory is 

= -V2^-^^^ , where = ^{x,k^) = ^ • 

The orbital angular momentum projection in this case is = —1. The spin structure 
indicated by perturbative theory provides a template for the numerator structure of 
the light-front wavefunctions even for composite systems. The structure of the elec- 
tron's Fock state in perturbative QED shows that it is natural to have a negative 
contribution from relative orbital angular momentum which balances the 5*2 of its 
photon constituents. We can also expect a significant orbital contribution to the pro- 
ton's Jz since gluons carry roughly half of the proton's momentum, thus providing 
insight into the "spin crisis" in QCD. 

The high x ^ 1 and high k± limits of the hadron wavefunctions control processes 
and reactions in which the hadron wavefunctions are highly stressed. Such configu- 
rations involve far-off-shell intermediate states and can be systematically treated in 
perturbation theory ||170| , |^. This leads to counting rule behavior for the quark and 
gluon distributions at x ^ 1. Notice that x —>■ 1 corresponds to k^ —oo for any 
constituent with nonzero mass or transverse momentum. 

The above discussion suggests that an approximate form for the hadron light-front 
wavefunctions could be constructed through variational principles and by minimizing 
the expectation value of H^c^ ■ 



12 Structure Functions are Not Parton Distribu- 
tions 

Ever since the earliest days of the parton model, it has been assumed that the leading- 
twist structure functions Fi{x,Q'^) measured in deep inelastic lepton scattering are 
determined by the probability distribution of quarks and gluons as determined by the 
light-front wavefunctions of the target. For example, the quark distribution is 

\lpn{Xi,k±i)\^'^6{xB -Xj). 

j=g 

The identification of structure functions with the square of light-front wavefunctions 
is usually made in LC gauge n ■ A = A'^ = 0, where the path-ordered exponential in 



Pq/N{XB,Q" 



Y[ dxi cPk^ 
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Figure 8: Glauber-Gribov shadowing involves interference between rescattering am- 
plitudes. 



the operator product for the forward virtual Compton amplitude apparently reduces 
to unity. Thus the deep inelastic lepton scattering cross section (DIS) appears to be 
fully determined by the probability distribution of partons in the target. However, 
Paul Hoyer, Nils Marchal, Stephane Peigne, Francesco Sannino, and I ||171|| have 
recently shown that the leading-twist contribution to DIS is affected by diffractive 
rescattering of a quark in the target, a coherent effect which is not included in the 
light-front wavef unctions, even in light-cone gauge. The distinction between structure 
functions and parton probabilities is already implied by the Glauber-Gribov picture 
of nuclear shadowing ||172| , |173| , |174| , p.75|] . In this framework shadowing arises from 



interference between complex rescattering amplitudes involving on-shell intermediate 
states, as in Fig. ||. In contrast, the wave function of a stable target is strictly real 
since it does not have on energy-shell configurations. A probabilistic interpretation 
of the DIS cross section is thus precluded. 

It is well-known that in Feynman and other covariant gauges one has to evaluate 
the corrections to the "handbag" diagram due to the final state interactions of the 
struck quark (the line carrying momentum pi in Fig. with the gauge field of 
the target. In light-cone gauge, this effect also involves rescattering of a spectator 
quark, the p2 line in Fig. ^. The light-cone gauge is singular - in particular, the gluon 
propagator d'^(~<{k) = -gr^^ —Q^" + "^^ntfc''"'' has a pole at /c"*" = which requires an 
analytic prescription. In final-state scattering involving on-shell intermediate states, 
the exchanged momentum is of O in the target rest frame, which enhances the 
second term in the propagator. This enhancement allows rescattering to contribute 
at leading twist even in LC gauge. 

The issues involving final state interactions even occur in the simple framework 
of abelian gauge theory with scalar quarks. Consider a frame with < 0. We 
can then distinguish FSI from ISI using LC time-ordered perturbation theory |^. 
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Figure 9: Two types of final state interactions, (a) Scattering of the antiquark (p2 
line), which in the aligned jet kinematics is part of the target dynamics, (b) Scattering 
of the current quark {pi line). For each LC time-ordered diagram, the potentially on- 
shell intermediate states - corresponding to the zeroes of the denominators Da, D^, Dc 
- are denoted by dashed lines. 



Figure ^ illustrates two LCPTH diagrams which contribute to the forward 7*T 
7*T amplitude, where the target T is taken to be a single quark. In the aligned 
jet kinematics the virtual photon fiuctuates into a qq pair with limited transverse 
momentum, and the (struck) quark takes nearly all the longitudinal momentum of 
the photon. The initial q and q momenta are denoted pi and P2 — ki, respectively. 

The calculation of the rescattering effect of DIS in Feynman and light-cone gauge 
through three loops is given in detail in Ref. ||171|| . The result can be resummed 
and is most easily expressed in eikonal form in terms of transverse distances r_i_,R± 
conjugate to p2±, k±. The deep inelastic cross section can be expressed as 



da 



where 



dQ^dxB 16772 y2 2Mv J P2 



\M{p^,f^,R, 



sm 



A{p2,f^,Rj 



g'W{f^,R^)/2 
is the resummed result. The Born amplitude is 

A{p2, f±, R±) = 2eg'MQp^ V{nnf^)Wir^, R^) 

where mf< = P2MXB + 'm? and 



(41) 



(42) 



(43) 



V{mrx) = 



1 



Ko(mr_L) 



(44) 



44 



The rescattering effect of the dipole of the qq is controlled by 



{2n) 



1 , \R±+r±\ 



(45) 



The fact that the coefficient of A in (|42|) is less than unity for all f±, R± shows that the 
rescattering corrections reduce the cross section. It is the analog of nuclear shadowing 
in our model. 

We have also found the same result for the deep inelastic cross sections in light- 
cone gauge. Three prescriptions for defining the propagator pole at /c"*" = have been 
used in the literature: 



1 

kf 



[kf - ir}i){kt + irji) 
kj - irji 
kf - ir]ie{k. 



-1 



(PV) 
(K) 
(ML) 



(46) 



the principal- value, Kovchegov ||1 76|| , and Mandelstam-Leibbrandt ||177|| prescriptions. 
The 'sign function' is denoted e(x) = G(x) — B(— x). With the PV prescription we 



have /„ 



Sdki 



0. Since an individual diagram may contain pole terms 



~ l/kf , its value can depend on the prescription used for light-cone gauge. However, 
the kf = poles cancel when all diagrams are added; the net is thus prescription- 
independent, and it agrees with the Feynman gauge result. It is interesting to note 
that the diagrams involving rescattering of the struck quark pi do not contribute to 
the leading-twist structure functions if we use the Kovchegov prescription to define 
the light-cone gauge. In other prescriptions for light-cone gauge the rescattering of 
the struck quark line pi leads to an infrared divergent phase factor exp i(f): 



9 



2 



Ko{XR±) + 0{g' 



(47) 



where A is an infrared regulator, and = 1 in the K prescription. The phase is 
exactly compensated by an equal and opposite phase from final-state interactions 
of line p2- This irrelevant change of phase can be understood by the fact that the 
different prescriptions are related by a residual gauge transformation proportional to 
6{k^) which leaves the light-cone gauge A'^ = condition unaffected. 

Diffractive contributions which leave the target intact thus contribute at leading 
twist to deep inelastic scattering. These contributions do not resolve the quark struc- 
ture of the target, and thus they are contributions to structure functions which are 
not parton probabilities. More generally, the rescattering contributions shadow and 
modify the observed inelastic contributions to DIS. 

Our analysis in the K prescription for light-cone gauge resembles the "covariant 
parton model" of Landshoff, Polkinghorne and Short [|1 78| , ^ when interpreted in 
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the target rest frame. In this description of small x DIS, the virtual photon with 
positive g"*" first splits into the pair pi and p2- The aligned quark pi has no final 
state interactions. However, the antiquark line p2 can interact in the target with an 
effective energy s oc k'j_/x while staying close to its mass shell. Thus at small x and 
large s, the antiquark p2 line can first multiple scatter in the target via pomeron and 
Reggeon exchange, and then it can finally scatter inelastically or be annihilated. The 
DIS cross section can thus be written as an integral of the a-qp^x cross section over the 
P2 virtuality. In this way, the shadowing of the antiquark in the nucleus a-qA^x cross 
section yields the nuclear shadowing of DIS ||174|| . Our analysis, when interpreted 



in frames with g+ > 0, also supports the color dipole description of deep inelastic 
lepton scattering at small x. Even in the case of the aligned jet configurations, one 
can understand DIS as due to the coherent color gauge interactions of the incoming 
quark-pair state of the photon interacting first coherently and finally incoherently in 
the target. 



13 A Light-Front Event Amplitude Generator 

The light-front formalism can be used as an "event amplitude generator" for high 
energy physics reactions where each particle's final state is completely labelled in 
momentum, helicity, and phase. The application of the light-front time evolution op- 
erator P~ to an initial state systematically generates the tree and virtual loop graphs 
of the T-matrix in light-front time-ordered perturbation theory in light-cone gauge. 
The loop integrals only involve integrations over the momenta of physical quanta 
and physical phase space '[\d'^k±idk^ . Renormalized amplitudes can be explicitly 
constructed by subtracting from the divergent loops amplitudes with nearly identi- 
cal integrands corresponding to the contribution of the relevant mass and coupling 
counter terms (the "alternating denominator method") [Q. The natural renormal- 
ization scheme to use for defining the coupling in the event amplitude generator is a 
physical effective charge such as the pinch scheme . The argument of the coupling 



is then unambiguous ||1 79|| . The DLCQ boundary conditions can be used to discretize 
the phase space and limit the number of contributing intermediate states without 
violating Lorentz invariance. Since one avoids dimensional regularization and non- 
physical ghost degrees of freedom, this method of generating events at the amplitude 
level could provide a simple but powerful tool for simulating events both in QCD and 
the Standard Model. 



14 The Light-Front Partition Function 

In the usual treatment of classical thermodynamics, one considers an ensemble of 
particles n = 1,2,...N which have energies {En} at a given "instant" time t. The 
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partition function is defined as Z = Z^nGxp— Similarly, in quantum mechanics, 
one defines a quantum-statistical partition function as Z = tr exp —j3H which sums 
over the exponentiated-weighted energy eigenvalues of the system. 

In the case of relativistic systems, it is natural to characterize the system at a 
given light-front time t = t + z/c; i.e., one determines the state of each particle in the 
ensemble as its encounters the light-front. Thus we can define a light-front partition 
function 



LC 



by summing over the particles' light-front energies p = — = xhe 



total momentum is P"*" = J2pt^ -^J- = J2nP-i-n, and the total mass is defined from 
P^P~ — P± = M^. The product -pLTic is boost invariant. In the center of mass 
frame where P = and thus P^ = P~ = M. It is also possible to consistently impose 
boundary conditions at fixed x~ = z — ct and x±, as in DLCQ. The momenta Pn,P±n 
then become discrete. The corresponding light-front quantum-statistical partition 
function is Z = tr exp —/SlcHlc where Hlc = is the light-front Hamiltonian. 

For non-relativistic systems the light-front partition function reduces to the stan- 
dard definition. However, the light-front partition function should be advantageous 
for analyzing relativistic systems such as heavy ion collisions, since, like true rapid- 
ity, y = \n^, light-front variables have simple behavior under Lorentz boosts. The 
light-front formalism also takes into account the point that a phase transition does 
not occur simultaneously in t, but propagates through the system with a finite wave 
velocity. 
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